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Abstract 

We compute the phase diagram in the N ^ oo hmit for lattice CP^^^ 
and QP^^^ ir-models with the quartic action, and more generally for mixed isovec- 
tor/isotensor models. We show that the N = oo limit exhibits phase transitions that 
are forbidden for any finite A^. We clarify the origin of these pathologies by examining 
the exact solution of the one-dimensional model: we find that there are complex zeros 
of the partition function that tend to the real axis as A^ ^ oo. We conjecture the cor- 
rect phase diagram for finite A^ as a function of the spatial dimension d. Along the way, 
we prove some new correlation inequalities for a class of A^-component cr-models, and 
we obtain some new results concerning the complex zeros of confluent hyper geometric 
functions. 
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1 Introduction 



Over the past decade there has been much interest in two-dimensional lattice cr-models 



taking values in the real projective space RP^ ^ 
I^, IT5|, the complex projective space CP^^^ p 



i, I, I, I, i i i I, H 
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12 




24 


or 
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quaternionic projective space QP^~^ (Hl-Q These models are good laboratories for investi- 
gating asymptotic freedom, topological phenomena and universality; they are intermediate in 
complexity between two-dimensional A^-vector models (cr-models taking values in the sphere 
5*^"^) and four-dimensional non-abelian gauge theories.^ 

In three dimensions, the RP^ lattice theory has long been employed as a model of nematic 
liquid crystals |3^, 0. In addition, a mixed isovector/isotensor model generalizing the 
RP^ theory has been proposed as a model of orientational transitions in crystals and 
of uniaxial liquid crystals exhibiting various kinds of head-tail asymmetry |36, 37, and 
a mixed isovector/isotensor model generalizing the RP^~^ theory arose in a study of N- 
component ferromagnets with annealed random couplings |^9[ 

The N ^ oo limit of the RP^~^ and CP^~^ cr-models has been studied by several groups 
in, 0, I, I, |, 0]; and Oku and Abe |3|, 0, Magnoli and Ravanini |[ and Ohno et al. 
have studied the N ^ oo limit of a mixed isovector/isotensor model that generalizes the 
jlpN-i modelj^ The result (in lattice dimension c? = 2) is that there is a first-order phase 
transition in the RP^~^ and CP^^^ models at /?c ~ 0.956. The question is whether the 
phase transition observed at N = oo occurs also at large but finite A^. 

The N = oo phase transition of the RP^~^ and CP^~^ models has two curious features 
0,0,1: 

(a) It occurs in all lattice dimensions, including d = 1. Indeed, it occurs even for a lattice 
consisting of only two sites! (However, in dimension d < 3/2 the transition is second- 
order rather than first-order f^.) 

(b) It is associated with the spontaneous breaking of the Z2 or f/(l) local gauge invariance. 
That is, in the mixed isovector/isotensor model for (3t > f^c, the free energy has unequal 
left and right derivatives with respect to Pv at /?y = 0. (Equivalently, the isovector 
energy Ey is a discontinuous function of /3v at = 0.) 

Both of these features are, however, rigorously forbidden at finite A^: in any short-range 
one-dimensional model (and of course on a finite lattice), the free energy and all expectation 
values are real-analytic functions of the parameters in the Hamiltonian^; and spontaneous 

^ There is also, of course, an enormous literature on the continuum versions of these models, which we 
refrain from citing. Our point of view is that the continuum model can be given a meaning only by defining 
it properly as a continuum limit of ultraviolet-cutoff models (such as lattice models) ; and in order to define 
such a continuum limit, one must first investigate the phase diagram of the cutoff models and search for 
critical points. 

^ These models are also good laboratories for the development and application of new collective-mode 
Monte Carlo algorithms, such as multigrid |^, and cluster algorithms ||, |3l[ |. 

^ We have also profited from consulting an unpublished manuscript on this problem that Sergio Caracciolo 
and Andrea Pelissetto have kindly shared with us [p^ . 

See El El and H, Theorem II. 5. 3 and Remark 1 following it]. 
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breaking of a local gauge symmetry cannot occur under any conditions [|^. Thus, it seems 
unlikely that such a transition can survive to finite A^. 

In this paper we propose to clarify this behavior, by comparing to the exact solution in 
the one-dimensional case.Q We shall see that, in d = 1, the phase transitions in the RP^~^, 
CP^-^ and QP^-^ models (and more generally in the mixed isovector/isotensor model) are 
artifacts of the N —>■ oo limit, arising from complex zeros of the partition function that tend 
to the real axis as ^ cxd. There is no phase transition for any finite We conjecture 
that the same is true in all lattice dimensions d <2. 

By contrast, in dimension d > 2 there is a transition from the isotropic phase to a phase 
with nematic long-range order this transition is expected to be first-order for 



all > 2 |5^. In this paper we shall compute the location and properties of this 

transition in the N ^ oo limit. 

The RP^~^, CP^~^ and QP'^~^ models to be considered here are defined by the lattice 
Hamiltonian 

where the sum runs over all nearest-neighbor pairs (xy) (each pair counted once), and each 
spin (Trc is a unit vector in M^, or Q^, respectively (in the real case, the complex 
conjugation and absolute value are of course redundant). The a priori measure for each spin 
cr^ is normalized uniform measure on the appropriate unit sphere, which we shall denote 
dQ{<Tx)- This model has a global symmetry group G = 0{N), U{N) or U{N, Q) as well as a 
local gauge group Gioc = Z2, U{1) or f/(l,Q) [~ SU{2)]. In order to treat all three cases in 
parallel, it is convenient to write K = M, C or Q, with corresponding real dimension k = 1, 
2 or 4. 

More generally, we shall consider the mixed isovector/isotensor model defined by the 
Hamiltonian 

_ ^2] 

2 



Pv (Tx ■ (Ty H —{(Tx-(Tyj 



:i.2) 



where each spin cr^. is a unit vector in M^. The pure RP^~^ model corresponds to /3y = 0. 
We shall restrict attention to the case Pt > 0. For simplicity, we refrain from considering in 
this paper the mixed models based on CP^~^ and QP'^"^. The central goal of this paper 
is to compute the phase diagram of the model ( |1.2| ) in the N ^ 00 limit — it turns out to 
be surprisingly intricate and dimension-dependent — and to discuss how much of this phase 
diagram should be expected to survive to finite N. 

^ The one-dimensional RP^~^ and mixed isovector/isotensor models have been studied previously by Liu 
and Joseph Bffl , Thorpe and BlumeB^ , Kohring and Shrock , Cucchieri et al. , Seller and Yildrim 



1)4 9[ , and Hasenbusch and Horgan |5C 

^ This phenomenon, in which a phase transition that is forbidden for all finite N nevertheless occurs 
at = 00, was observed nearly two decades ago by Celmaster and Green jsi). They noted that the one- 
dimensional SU (N) and U (N) principal chiral models at A'^ = cxd show spontaneous breaking of a continuous 
global symmetry, and that the two-dimensional SU{N) and U{N) lattice gauge theories at = cx3 show 
spontaneous breaking of the local gauge symmetry — transitions that are of course forbidden at finite N. 
They concluded that these examples "illustrate the possible pitfalls of investigating symmetry breaking at 
N 00." 
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The plan of this paper is as follows: In Section ^ we show how to obtain the N oo limit 
(and the expansion) for the mixed isovector/isotensor model ( |1.2|) , paying particular 
attention to conceptual questions (the order of integration, how to determine which saddle 
point is dominant, ...). In Section ^ we specialize to the translation-invariant nearest- 
neighbor case, and give two complementary methods for computing the N ^ oo limit: a 
direct method and a parametric method. In Sections ^ and ^ we carry out these calculations 
in detail for /3y = and Pv 0; respectively; our main goal is to determine the phase diagram 
in the {jSy, /3T)-plane as a function of the spatial dimension d (which we treat as a continuous 
variable). Some of the results in these sections are known from previous work 0, ^, but 
most of them are new. In Sections |^ and we solve the one-dimensional model and extract 
the N oo limit; we clarify the origin of the phase transition at = oo, and explain why it 
does not occur for any finite N. In Section we draw some conclusions, and conjecture the 
phase diagram in the (/3y, /3T)-plane for finite A^. In Appendix ^ we collect some formulae 
concerning the lattice propagator and free energy in general dimension d. In Appendix P we 
prove some asymptotic expansions for confluent hypergeometric functions and their zeros; 
these appear to be new. In Appendix |C| we prove some new correlation inequalities for a 
class of A^-component a- models that includes the mixed isovector/isotensor model. 



2 N ^ oo Limit for Mixed Isovector/isotensor Models: 
General Theory 



Let us begin by considering the mixed isovector/isotensor model on an arbitrary finite 
lattice C having V sites. The Hamiltonian is 



H 



K 



2 ■ ^y) 



(2.1) 



where the sum runs over all unordered pairs (xy) of sites in C (each pair counted once). We 
shall assume that K^y = Kyx > 0, while J^y = Jyx may have either sign. Without loss of 
generality we assume a priori measure for each spin cr^ is normalized 

uniform measure on the unit sphere in R^: 



dVL{(T^) 

The partition function is therefore 



r(A^/2) 



5{(Tl-l)d^(T^. 



TT 



N/2 



(2.2) 



r(iv/2) \ 




Jxy^x ' H~ 



K. 



xy 



(^x-(Ty) 



\[5{cTl-l)d^a,. (2.3) 



Our main goal is to compute the leading term in the large- A^ expansion of the free energy 
— log Z, namely the limiting free energy per component 



— lim — log Z . 



(2.4) 
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We also want to compute the leading term in the large- expansion of the n-point correla- 



tion functions (cxi"^^ ■ 
correlation functions 



cxi""'*), the most important of which are the isovector and isotensor 



Gv{x,y) 



(2.5) 



N 



GT{x,y) = (T.-T,) ^ ^ (Ti"^)T(-^)) 

a,l3=l 
1 



y> , ^ 

where we have defined the symmetric traceless isotensor field 



(2.6) 



(2.7) 



We begin by introducing auxiliary variables iixy so as to "Gaussianize" the isotensor part 
of the Boltzmann weight: 



exp 



xy 



-1/2 



exp 



l^xy 



2K 



+ fJ'xyf^. 



X '-' y 



xy 



xy ■ 



(2.8) 



xyj 



-1/2 



exp[—fil /2Kxy] is of course interpreted as 5(/i, 



xyJ 



[When Kxy = 0, the Gaussian {27cK. 
It follows that the partition function ( p.3|) of the mixed isovector /isotensor model can be 
rewritten as that of a pure A^-vector model with annealed random couplings J^y + fixy, where 
the fixy are Gaussian-distributed with mean zero and variance K^yfl 

Next we introduce auxiliary variables {a^} in the standard way to "Gaussianize" the 
a priori measure: 



' 2m J 



(2.9) 



where = Cx{{fi}) is a real constant chosen large enough to make the integral over {cr} 
absolutely convergent. We then get 



exp 



f^xy 

2K. 



xy 



where 



B 



xy 



20ix^xy i.'^xy ~l~ f^xy) • 



(2.10) 



(2.11) 



^ This equivalence was noted almost two decades ago by Oku and Abe ||3^ , who started from the A'^- vector 
model with annealed random couplings and showed its equivalence with the mixed isovector /isotensor model 
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Note that the order of integration here is crucial, as the {cx} must be chosen in a {/i}- 
dependent way so that Re 5 is positive definite. 

We can now perform the Gaussian integral over {cr}. For convenience let us introduce new 
variables J^y = Jxy/N, K^y = K^y/N, Jl^y = fixy/N, = a^/N and hence Bxy = Bxy/N. 
We get 

(2.12) 

where 

T{^,a) = -^5,. + J]-^^ + itrlog5 + ^. (2.13) 

X (xy) ^^xy 

Let US now perform the {5} integrals by the saddle-point method. The saddle-point 
equations are 

^ = {B-^).. -1 = for all X . (2.14) 
It can be proven that there exists a unique real saddle point 5^ = 5^ (/I) in the region 



where B is positive-definite; we conjecture (but have not proven) that this is the only saddle 
point, real or complex, in the region where ReS is positive-definite. For the purposes of 
calculating the leading — oo contribution, it suffices to replace 5 by its saddle-point 



= J^{Jl,a,{Jl)) . (2.15) 

We then perform the integral over {Jl} by Laplace's method: it is dominated by the config- 
uration(s) of {Jl} that make J-'{J1) an absolute minimum, and so we have 



inf T{J1) (2.16) 



(the prefactors cancel to leading order). 

The stationary points of J^{J1) are given by 



& - {B-%y = for alia;,!/. (2.17) 
K 



In general, !F has many stationary points. Indeed, for a pure RP^~^ model (J = 0), any 
gauge transform Jlxy VxVyJ^xy iv = ±1) of a stationary point is also a stationary point 
with the same value of JF, so in particular there are many degenerate absolute minima. If 
Jxy ^ 0, there is no longer a gauge invariance, but at least for small Jxy all of the foregoing 
stationary points survive (and move slightly). We conjecture that if Jxy > for all x,y 
("ferromagnetism"), then one of the absolute minima has Jlxy > for all x, y, and moreover 

^ In principle one should prove that the limit N ^ oo and the integration over {/!} can be interchanged. 
Unfortunately, we have no rigorous proof, though it seems eminently reasonable. 
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there is a unique absolute minimum with this property; but we have, as yet, no proof. Let 
us also remark that if J = 0, then /I = is always a stationary point (with 5*^. = | for all 
x), but this is the dominant saddle point only for K sufficiently small (see Section 

By expanding ( |2.12| ) / ( ^.13| ) around the saddle point, one can obtain a systematic expan- 
sion in powers of Examination of this expansion could potentially shed alternative 
light on the reliability of the large- limit; but we shall not pursue this approach here. 

The saddle-point method also provides a 1/N expansion of the correlation functions. At 
leading order in the fields cr = {(ji°'^}xec,i<a<N are Gaussian with mean zero and 

covariance matrix 

= ^ [iB;%y + 0(l/iV)] (2.18) 

where 

^i^xy ^O^^x^xy ^Jxy ~t~ f^xy^ • (2.19) 

Thus, 

Gv{x,y) = {(Tx-(Ty) = iB-')xy + 0{l/N) , (2.20) 

where the saddle-point equation (|2.14|) ensures consistency when x = y (since cTj, • cTj, = 1). 
Moreover, higher-point correlation functions are given by sums of products of B~^: for 
example. 



{a2) (as) (04) 



^t^^^x^^) = 175 {B-^)x,x,{B-^)x,x, + two permutations + 0{l/N) 



'xi " X2 " X3, '^X4 I 

(2.21) 

and likewise for all 2n-point correlation functions. One useful consequence is that, to leading 
order in 1/A^, correlations involving products of 0(A^)-invariant combinations factorize 



n '^-^ ^^y ) = ^""-^ ■''y^) + 0{i/N) . (2.22) 



ij=l / i=l 

In particular, 

GT{x,y) ^ {{cTx ■ cTyf) - 1 = {B:')ly + 0(1/Ar) . (2.23) 

For the pure RP^~^ model, the factorization ( p.22| ) should be understood as applying to the 
symmetry-broken measure {■)sb obtained by letting /5y | 0. Alternatively, for Z2-gauge- 
invariant correlations — that is, those for which each lattice site occurs an even number of 
times in the sequence xi, . . . , Xn,yi, ■ ■ ■ ,yn — we have (to leading order in 1/A^) 



\2\l/2 

where ( ■ ) denotes expectation in the Z'2-gauge-invariant Gibbs measure. 
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N ^ oo Limit for Mixed Isovector/Isotensor Models: 
Translation-Invariant Case 



3.1 Basic Definitions 

Let us now restrict attention to the case where £ is a d-dimensional cube of side L 
with periodic boundary conditions, and the couphngs are translation-invariant and nearest- 
neighbor: J^y = PvS^y, Kxy = Prf^^y with Pt > 0. Let us also restrict attention to 
solutions for {5} and {]!} that are translation-invariant and for which Ji is the same in all 
lattice directions. (In principle we should prove that these are the dominant saddle points; 
but since we are unable to do so, we shall simply assume that they are, like all previous 
workers!) We thus have 

B,y = 2aS,y - 0v + J1)S!''' 



or in Fourier space 



where 



Bip) 



25 



7^ COS Pi 

i=l 



7 

Let us introduce the fundamental quantities 



a 



9d,L{l) 



1-7E 



COS Pi 



i=l 



COS Pi 



(3.1) 
(3.2) 

(3.3) 

(3.4) 
(3.5) 



where the sum runs over momenta with components pi = 2Tiki/ L {ki = 0, 1, ... ,L — 1). 
These quantities are well-defined for |7| < 1/d; they are strictly increasing functions of |7|, 
and diverge as |7| t 1/d. They satisfy 



We then have 



d9d,L{l) 



^Al(7) 



7 



1 



2(i«7 



[fd,M - 1] for \x-y\^l 



(3.6) 

(3.7a) 
(3.7b) 



where of course the right-hand side of (3.7b) must be interpreted as zero when 7 = 0. After 
passing to the infinite-volume limit, the correlation function {B^^)xy exhibits exponential 
decay 



xy 



~ e 



-m\x-y\ 



as \x — y\ 



oo 



(3.8) 
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with mass gap 



--{d-1) 

,7 



m = arccosh 

Our goal is to compute the free energy per unit volume 

F = L-^F 



(3.9) 



(3.10) 



in the infinite-volume limit L —>■ oo. We shall do this in two ways: a direct approach (Section 



3^ and a parametric approach (Section p.3|) . These two approaches have complementary 
advantages, and together they give added insight into the phase diagram. The parametric 
approach has been used previously by Magnoli and Ravanini 0, while the direct approach 
appears to be new. 



3.2 Direct Approach 

We implement the method set forth in Section |^: namely, we compute the function 
J-'iJl) = J-'ip, 5^(/^)) and then find its absolute minimum. The first step is to fix Jl and solve 
the saddle-point equation (|2.14| ) for 5: from ( |3.7a| ) we have 



= fdAl*) = fd,L\ 



(3.11) 



It is not hard to see that 7^ goes monotonically from —1/d to 1/d as Jl goes from —00 to 
00 .0 So we can use 7^ (henceforth denoted simply by 7) as an independent variable in place 
of yU. The free energy per unit volume J-'iJl) = L~'^J^(J1), written as a function of 7, is 



1 - fdAi) + 



(3i 



- pv] + log fdAi) + 9dAA 



Let us write it as 



d 



HA = -^0(7) + /?y^i(7) + 
where ^^0(7) [resp. ^^1(7)] is an even (resp. odd) function of 7: 

d 



V 1 



^0(7) 



-^1(7) 



1 - fdAA + -^AfdAA^ + logfdAA + 9dAA 

4(3t 



d 



2Pj 



ifdAA 



The free energy per unit volume is then 



F = inf^(7) 

7 



(3.12) 
(3.13) 

(3.14) 
(3.15) 

(3.16) 



^ This is easiest seen in reverse: write /x = \'^*fd,L{"f*) — Pv- It follows from the properties of fd,L that 
Jl goes monotonically from —00 to 00 as 7* goes from —1/d to 1/d. 
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Since J^i(7)/7 < for all 7, it follows that when /5y > (resp. /3y < 0) the absolute minimum 
of T lies at 7 > (resp. 7 < 0). Of course, when /?y = 0, is an even function of 7, so the 
sign of 7 is irrelevant. 

Note also that fd,L{'j) goes monotonically from 1 to +00 as I7I goes from to l/rf. So 
we can use / = fd,L{l) and sgn7 as independent variables in place of 7, if we prefer. As we 
shall see, this is sometimes convenient. 

Note, finally, that we can read off from the /3y = solution the qualitative behavior of 
the free energy for small but nonzero (3v- Indeed, suppose that the absolute minima of ^^0(7) 
lie at ±7i, . . . , ±7n, with common minimum value jFo(±7j) = Fq. Then, for small (3vi we 
have from (|3.13|) 



F = inf ^(7) = Fo - \Pv\ max |J-i(7,)| + OiP^y] 

7 l<j<n 



(3.17) 



Thus, if the unique absolute minimum of ^^0(7) lies at 7 = 0, we have F = Fq + 0(/3y); but 
if there is any absolute minimum of ^^0(7) at 7 7^ 0, we have F = Fq — 0{\Pv\) and hence a 
first-order phase transition as jSy passes through 0. Surprisingly, this happens even in finite 
volume (L < 00). Note also that the isovector energy Ey ~ dF/dPv is discontinuous at 
such a first-order transition, indicating that the Z2 gauge invariance of the jSy = theory is 
spontaneously broken. Both these features are manifestly artifacts of the N ^ 00 limit, as 
they are rigorously forbidden at any finite N. 

Let us now consider the infinite-volume limit L 
infinite- volume limit in (|3.4]) and (p.5|), yielding 



00. Formally we can pass to the 



fdil) = lim fdAl) 



9d{l) = lim gd,L{l) 



cos Pi 



log 



7^ COS Pi 

i=l 



(3.18) 



(3.19) 



where the integral runs over the Brillouin zone [— vr,??]'^. This substitution gives the correct 
answer in dimension d < 2. However, in dimension d > 2 we must be a bit more careful, 
because the infinite- volume quantity fd{,l) tends to a finite value as 7 j l/d even though 
the finite- volume quantity fd,L{l) diverges due to the contribution of the p = mode.|3 In 
particular, the equation fd,L{l) = f has a solution for arbitrarily large /, which satisfies 
7 = 1/d — 0{L~'^). In infinite volume, we can think of all these solutions with / > /* as 
arising from 7 = l/rf; we shall refer to these solutions as the "condensate phase" , and use / 
to parametrize them rather than 7. It is important to note that when / = 25^ exceeds f^, 
the finite-lattice Green's function B^{p)~^ defined by (|3.2|) , which consists of delta functions 
located at the momenta with components pj = 27rA;j/L, converges in distribution as L ^ 00 
not to the corresponding infinite-lattice Green's function B^{p)~^, but rather to 



k 

f 



6{p) 



(3.20) 
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The same happens as 7 | — l/d, due to the contribution of the p = (tt, . . . , tt) mode. 
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Here the delta- function contribution at p = corresponds to the long-range order 



hm {b:': 



\x-y\^oo 



xy 



-, /* 



(3.21) 



Note also that the zero mode does not contribute to gd,L in the infinite-volume limit at fixed 
/ > /*, since l/rf— I7I > L"'^ and limi^oo(logi^'^)/i^'^ = 0; therefore we should consider 
9d,L{l) to be fixed at the value = gd{l/d). 

In summary, in dimension d > 2 we should consider the additional branches 

^±(/) = ^o*(/) ± Pv^Fuif) + 



2A 



(corresponding to 7 = ±l/d) as a function of the variable / G [/*, 00), where 



^o*(/) 



f + 



+ log/ + g^ 



(3.22) 



(3.23) 



(3.24) 



and we should then minimize over both (|3.13| ) and ( p.22| ). Since J-'u{f) < 0, the optimal 
choice of sign in (|3.22| ) is clearly sgn7 = sgnf3v Whenever the absolute minimum of the 
free energy is found in the condensate phase / > /* at 7 = l/rf, the correlation functions at 
N = 00 exhibit both ferromagnetic and nematic long-range order^: 



lim Gvix,y) = ^ 

\x—y\^oo J 



lim Grix^y) 



"7 



(3.25) 



(3.26) 



In Appendix ^ we collect some results concerning the behavior of the functions /d(7) 
and gd{l)- In particular, we show how these quantities have a natural analytic continuation 
in the parameter which allows us to discuss noninteger dimensions ci > 0; the ability to 
treat ci as a continuous variable is very useful. The small-7 Taylor expansions of fd{l) and 
5^(7) are given in (|A.13|) and (|A.14| ), respectively. It follows that ^0(7) and ^1(7) have the 
following small-7 expansions: 



•^0(7) 



d 



'T 



1 y + 



d fM 



64 



+ (3 - lOrf) 7 



+ 



d f6d{7d-3) 



192 V 



{62d^ - 63d + 20) 7^ + 



^1(7) 



d 



d^ 



-7 



-7 



3d\2d 



'-7' + 



2Pt 4:pT 16/5t 
The behavior of /d(7) and gdil) as 7 t l/rf is discussed at the end of Appendix |A.l 



(3.27) 
(3.28) 



In the antiferromagnetic case 7 — — 1/c?, one has instead hm (—1)'^ ^'G'y(x,j/) = 1 — /*//. 

\x-y\^oo 
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3.3 Parametric Approach 



Thus far we have been analyzing the phase diagram by searching directly for the absolute 
minimum of J-'iJi), without bothering to use the stationarity equation (jl) = 0. On the 
other hand, in searching for the absolute minimum we have the right to restrict attention 
to such stationary points if we so choose. As we shall see, this alternate approach has some 
advantages: in particular, it will allow us to visualize the phase diagram simultaneously for 
all /5t > (at any fixed value of Pv)- 

We proceed as follows: The stationarity equation ( |2.17|) for /I, when combined with 
(p.3|)/(3.7b) and the saddle-point equation ( |3.11| ) for 5, yieldsQ 



(3.29) 



If we restrict attention to solutions of ( p.29|) , then the free energy (|3.12| ) can be rewritten so 
as to eliminate reference to either /3y or I3t'- 



1 - fdAi) + iog/rf,L(7) + SdAi) + 



/5t //d,L(7)-l 



d \ ifdAl) 



^ + log/d,L(7) + 9d,L{V - Pv 



ifdAl) 



(3.30a) 



(3.30b) 



The solutions to (|3.29| ) come in two types: 

(a) 7 = is a solution to ( p. 291) if and only if /Jy = (i.e. the pure RP^~^ model); the 
corresponding free energy is JF = 0.Q ^ 

(b) If 7 7^ 0, we can solve (jpl) for /3t: 



dl'fdAlf 



d^fdAl) 



V ■ 



(3.31) 



For each fixed equations ( p.31|) and (3.30b) give [3t and parametrically as a function 
of 7 G (— l/(i, 1/d). We must trace out this curve in the (/5t, JF)-plane, and then take the 
smallest value of for each /?t-|3 

In Sections ^ and |] we will carry out this procedure (in the limit L oo) and discuss the 
qualitative properties of the solution as a function of the lattice dimension d. As preparation. 



Alternatively, this same equation can be obtained by differentiating (3.13) with respect to 7, using (p^q), 
and setting the result equal to zero. Miraculously, the derivative of fd,L enters here only through the overall 
factor 7/(j,L + f'd l [which never vanishes, see e.g. ( A.IOD ], so the solution does not involve ^ at all. 

13 remark that Magnoli and Ravanini equation (3)] made the mistake of multiplying equation ( 3.29| ) 
by 7, thereby introducing a spurious 7 = solution also for (3v 7^ 0. Fortunately, this error had no effect on 
their final results, because the purported 7 = solution was never the dominant saddle point anyway. 



^'^ Alternatively we can fix /3t and then use equations (3.31) and (3.30a) to give /Sy and parametrically 
as a function of 7. Yet another variant was employed by Magnoli and Ravanini who fixed p = Pt/Pv 
and computed f3 = (iy + (3t and J- parametrically as a function of 7. 
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let us compute, using / ( [A78|) , the derivatives with respect to 7 of (3.30b) and (|3.31| ): 

(3.32) 



7(/-2)/' + 2/(/-l) ^ f-f + ^f 
Pv 



47/ 

6;7/[7(/-2)r + 2/(/-l)] 
2(/-l)2 

rf7(2^y-27/ + 7/'; 



+ Pi 



272/2 
d{f -P + if) 



(3.33) 



2(/-l)2 



r + 



if -If ^ 



^ 2d(/3y-27/ + 7/^) ^; I dp 



/-I 



(3.34) 



where /' = df /dj; these equations hold both in finite volume and in the infinite- volume 
limit. In particular, we have 



dj" 

d^T 



2d^^P 



< . 



(3.35) 



It is a curious fact that the derivative dJ^/d^T does not depend on jSy- Note also that 
dT jd^T becomes more negative I7I grows (in all dimensions d > 0): this strict monotonicity 
follows from inequality ( [A.10[ ). 

In the infinite- volume theory in dimension c? > 2, we need to consider also the "condensate 
phase" parametrized by / G [f^, 00) with 7 = (l/rf) sgn Pv and g = Here (3.30) is replaced 
by 



1- f + log/ + g, + dPr 
1 - f 

— ^ +logf + - d\pv\ 



/-I 
/ 

/-I 
/ 



while ( |3.31| ) is replaced by 



Pi 



P 



f 



2d{f-l) /-I 
The derivatives of (3.36b) and ( |3.37| ) with respect to / are 



\Pv\ . 



and hence 



dT_ 
W 

dPr 
W 

dT 

dPr 



p-2f + 2d\P 



v\ 



4/2 

p-2f + 2d\P, 
2d{f - 1)2 

d{f-V^ 



2f 



< 



(3.36a) 
(3.36b) 

(3.37) 

(3.38) 
(3.39) 

(3.40) 
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Once again, dJ^/d[3T is independent of (3v- Moreover, (|3.35| ) and (|3.40|) coincide in the 
condensate phase I7I = l/c?.Q In particular, they coincide at their joining point I7I = l/d, 
/ = this means that any phase transition between these two branches must be of second 
(or higher) order. Note also that dJ^/d(3T is a monotonic function of /: it becomes more 
negative as / grows. 

For (i > 2 it is convenient to use / as the independent variable in both the "normal 
phase" (I7I < and the condensate phase: see e.g. Figure |^ below. In particular, 

equations ( |3.32| )/ (|3.33| ) for the normal phase can be converted to yield dTjdj and d(3T/df 
by dividing them by /'(7). It is curious to note that in dimensions 2 < < 4, where 
fl{d) = lim dfd{'y)/d'y is infinite, the limiting slopes dj^/df and d^r/df for 7 | l/ci in the 

normal phase coincide with the values (|3.38|) / (p.39|) obtained for / | in the condensate 
phase (this matching of slopes is seen e.g. in Figure |^). On the other hand, for ci > 4 this 
matching of slopes does not occur. 



4 Behavior of Solutions for = ^ 
4.1 Overview 

Our conclusions concerning the phase diagram of the N = 00 RP^^^ model, as a function 
of the spatial dimension d (treated as a real variable), are summarized in Figure |1|. There are 
three possible phases: a white-noise phase (7 = 0) in which the isotensor correlation length is 
zero; an isotropic phase (0 < I7I < 1/d) in which the isotensor correlation length is finite but 
nonzero; and a nematic ordered phase (I7I = 1/d) in which the isotensor correlation function 
exhibits long-range order. For < d < 3/2, there is a second-order transition at /?r = 1 from 
the white-noise phase to the isotropic phase. For 3/2 < d < d^ ^ 2.38403, this transition 
is first-order and takes place at a dimension-dependent (3t,c < 1- For 2 < d < d^, there is a 
subsequent second-order J;ransition from the isotropic phase to a nematic ordered phase at 
a dimension-dependent Pt,c' = f^/['^d{f* — 1)] < 00. At d = d^ the two transitions collide 
at /3t,* ~ 0.873264, and for d> d^ there is only a first-order transition from the white-noise 
phase directly into the nematic ordered phase, at a dimension-dependent /3t,c" < 1-0 

Let us now show how these conclusions can be derived using either the direct or the 
parametric approach. 



4.2 Direct Approach 

For [3y = 0, the free energy is given simply by J-'oij)- Since JFq is an even function of 7, 
one of its stationary points is always 7 = 0. To learn about other possible stationary points, 
we begin by examining the Taylor expansion ( |3.27D of ^^0(7) near 7 = 0. The coefficient 



Both these properties foUow alternatively from the fact that ( |3.35 ) holds also in finite volume 



Ohno et al. ^ were the first to note the existence of the special dimension d.^, ; but they did not attempt 
to determine it numerically, stating only that 2 < d^,, < 3. Caracciolo and Pelissetto, in an unpublished 
manuscript [n2[, found essentially the same value of as we report here. 
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of 7^ is positive for /3r < 1 and negative for /5t > 1- The coefficient of 7^ is positive for 
(3t < 8d/ max(10(i — 3, 0) and negative for (3^ > 8d/ max(10(i — 3, 0). It follows that the local 
behavior near 7 = falls into three cases: 

(a) < d < 3/2: The coefficient of •y^ is positive when /3r = 1- Thus, as Pt passes through 
1, the local minimum at 7 = turns into a local maximum, and a pair of degenerate 
minima at 7 = ±7o(/3t) bifurcates off it. These minima satisfy the equation 



/3i 



2[Wo) - 1] • 



(4.1) 



As /3t i 1 we have 



T 



l-d 



1 + 



d 



1/2 



a/2 



(Pt - 1, 
20c/2 - 54rf + 31 



X 



Ail-df 



{Pt-1) + OiiPr-l) 



•T 



8(1 -d) 

+o(0T-ir 



rf(4c/2 _ I8rf + 13) ~ 3 
48(1 -ci)3 



(4.2) 



(4.3) 



(b) d = 3/2: The coefficient of 7^ is zero when /3t = 1, but the coefficient of 7^ is positive. 
Therefore, the qualitative behavior is the same as for d < 3/2, but the powers are 
changed: 

(4.4) 

:Fo{1o0t)Jt) = -^0T - if - - ir + O[0T - If) 

(4.5) 

(c) d > 3/2: The coefficient of 7^ is negative when /5t = 1. The behavior for /5t ~ 1 
and 7^0 will turn out to be irrelevant, as a first-order phase transition occurs before 
/3t = 1; at a 7 far away from 0. 

The local behavior near 7 = does not, of course, necessarily determine the global 
behavior of ^o(7)-Q But our numerical computations confirm — although we have been 
unable to prove it analytically — that ^^0(7) always has at most one pair of local minima 
7 = ±7o(/3r) 7^ 0, and that: 

This point was overlooked by Hikami and Maskawa ^ Appendix], who asserted (incorrectly) that there 
is a phase transition at fix — I in all dimensions d. 
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(a) < (i < 3/2: For /3t < 1, the unique absolute minimum of J^o(7) lies at 7 = 0. For 
/3t > 1, ^0(7) has precisely two degenerate absolute minima, at 7 = ±70 (/?t) satisfying 
( [4 .21 ). The free energy F therefore equals for /5t < 1, and it satisfies ( [4 .31) for /3t > 1. 
There is thus a second-order phase transition at /5r = 1, in which the specific heat 
Ch ~ d'^F/dP^ has a jump discontinuity (see Figure ^). 

(b) d = 3/2: The quahtative behavior is the same as for d < 3/2, but 7o(/5r) satisfies 
( [1.4|) , and the free energy for /3r > 1 satisfies (|4.5|) . There is thus a second-order phase 
transition at /3t = 1, in which the specific heat has a (/3r — 1)"^'^^ singularity as /3t i 1 
(see Figure |]b). 

(c) 3/2 < (i < 2: For /5t < a certain /3r,c(c?) < 1, the unique absolute minimum of ^0(7) 
lies at 7 = 0. For Pt > PtA^)^ ^oil) has precisely two degenerate absolute minima, 
at 7 = ±Jo{Pt), where 



7c > as Pt i Pt,c 
1/d as Pt ^ +00 



(4.6) 



There is thus a first-order phase transition at a dimension-dependent PT,c{d) < 1 (see 
Figure ||c for the case ci = 2).0 We have 



Ic 



1 - A(d 



3n2 
2> 



100 



(d- 



o((rf-r 



(rf- 1)1/2 + |(rf- 1)3/2 + 0({d- 



3N5/2 



(4.7) 



(4i 



as i 3/2. 

For (i > 2 we need to consider also the "condensate phase" (|3.23|) parametrized by / > 
at 7 = Ijd and g = g*. For convenience we use / as the independent variable in the plots 
also in the "ordinary phase" < 7 < l/rf. We have: 

(d) 2 < (i < (i* ~ 2.38403: The qualitative behavior near the first-order transition point 
Pt,c is the same as for d < 2. However, beyond this point the behavior is different: 
loiPr) now reaches 1/d when Pt reaches the finite value 



PT,Ad) = 



2d{f, 



(4.9) 



For Pt > Pt,c', the function ^^0(7) is strictly decreasing on the interval [0, l/d], and 
the absolute minimum of the free energy is found in the condensate phase at 

1/2- 

1+11-^1 > /, 



/ 



dPT 

/-I 

4 



+ 



dpT 
log f + g* 



< 



For d = 2, the transition at /3t,c ~ 0.956 was found previously by Magnoli and Ravanini 



(4.10) 

(4.11) 

using the 



parametric approach (see Section 4.3 below). 
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There is a second-order isotropic-to-nematic phase transition at [3t = Pt,c' (see Figure 
^). The two transitions Pt,c < Pt,c' come closer together as d grows, and they colhde 
at the dimension where 7c reaches 1/d (see Figure^). Equivalently, is defined by 
the condition 

r 1 1 . , r < for 2 < c/ < rf* 
= -^ + i^^^4^ =0 for d = d. (4.12) 



> for d> d^ 



(see Figure ^. 



(e) d = d^: The two transitions coincide: /?t,c = Pt,c' = Pt,* ~ 0.873264. Thus, for 
Pt < Pt,*, the unique minimum hes at 7 = 0; for Pt = Pt,*, it lies exactly at 7 = l/ci 
with f = f*] for /5t > it lies in the condensate phase 7 = l/rf at the value ( |4.1CI| ). 
There is a first-order phase transition as Pt passes through Pt,* (see Figure ^). 

(f) d > d*: For Pt < f^/[2d{f* — 1)], J^o(7) is a strictly increasing function of 7 on [0, l/d]; 
for f^/[2d{f* — 1)] < /^■r < 1, it first increases and then decreases; for Pt > 1, it is 
strictly decreasing. In all three cases, ^0(7) never has a local minimum in the interval 
< 7 < l/(i. Therefore, the absolute minimum is either at 7 = or else in the 
condensate phase at ( [4.10|) /( ^.11|) . There is a PT,c"{d) such that the former behavior 



occurs for Pt < PT,c"{d) and the latter for Pt > PT,c"{d) (see Figure ^jc): it is given by 

where fc" is the unique solution of log / + (?* = (/ — l)/2 satisfying / > As 
c? — > 00, we have fc" = foo + 0{d^^) and Pt,c" = /^/[2(/oo — 1)]*^""^ + O^d^"^), where 
/oo ~ 3.512863 is the unique solution of log/ = (/ — l)/2 on / > 1. 

4.3 Parametric Approach 

Since and Pt are even functions of 7 when Pv = 0, it suffices to consider 7 > 0. From 
( p.31| ) we have 



~ d^^P 1^3-2^2 20^^-54^ + 31 4 , 

Pt = = 1 + -^1 7 + ... . (4.14) 

In Figure ^ we plot Pt versus 7 for dimensions (i = 1, 3/2, 2, 9/4, ci^,, 3. Moreover, for c? > 2 
we need to consider the "condensate phase" 

parametrized by / G [/*, 00). 

The following assertions appear to be true: 

Since — (/— l)/4+ (log/ + (7*)/2 is a concave function of / that is > at / = /* and tends to —00 as 
/ +00, it has a unique zero in the interval (/*, 00). 
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(a) For < (i < 3/2, is a strictly increasing function of 7 on [0, 1/d). It tends to +00 
as 7 I 1/d. 

(b) For 3/2 < d < d ^ 2.55391, there exists •jcuspid) G (0, 1/d) such that Pt is a strictly de- 
creasing function of 7 on [0, •jcusp] and a strictly increasing function of 7 on [7cnsp, 1/d). 
The value of 'jcusp is determined by the condition 7 = 2/(/ — l)/[(2 — /)/']. Moreover, 
d'Jcusp is an increasing function of d: it moves from jcusp = 0at(i = 3/2to 'jcusp = 1/d 
at d = d. The value of d is determined by the condition = 2. 

(bi) For 3/2 < < 2, tends to +00 as 7 t l/d. 

(b2) For 2 < d < d, Pt tends to a finite value f^/[2d{f^ — 1)] as 7 t l/c^- However, 

the slope d^r/dj at this point is +00. In the condensate phase / > Pt is a 
strictly increasing function of /. 

(c) For d > d, Pt is a strictly decreasing function of 7 on [0, 1/d). It tends to a finite value 
/2/[2f/(/* - 1)] < 1 as 7 T l/d. 

(ci) For d < d < 4, the slope dPr/d'j is —00 at 7 = l/ci. 
(C2) For (i > 4, the slope is a finite negative number. 

In either case, in the condensate phase / > /*, at first decreases and then increases; 
it reaches its absolute minimum at / = 2 and Pt = 2/d. 

The assertions about the behavior as 7 t 1/d are simple rigorous consequences of known 
facts about fdi'j) [see Appendix]^, and the assertions about the condensate phase are trivial 
calculus. But we have been unable to demonstrate analytically the observed monotonicity 
properties of Pt{i)- 

Let us next look at the free energy ^-"(7) defined by (3.30b) with Pv = 0: 

_ log/ + ^ ^ d{2d~3) 4 d{7d^-18d+10) 6 

4 2 64 ^ 48 7 ••• • I • j 



From (|3.35 ) it follows that JF is decreasing wherever Pt is increasing, and conversely. As 



7 I l/d, JF tends to —00 [resp. to a finite number = — (/* — l)/4 + (log/* + g*)/2] in 
dimension d <2 (resp. d > 2). In dimension c? > 2, we need to consider also the condensate 
phase / > /*, in which JF is given by (3.36b). Again we have to distinguish 2 < d < d 
(corresponding to /* > 2) from d > d (corresponding to < 2). In the former case, J-' 
decreases monotonically from JF* to —00 as / rises from to +00. In the latter case, JF first 
rises from JF* to (2 log 2 — l)/4 + g^/2 as / rises from to 2, then decreases monotonically 
to —00 as / rises from 2 to +00. We refrain from showing the plots of JF versus 7. 



We are now ready to plot JF parametrically versus Pt- As in Section |4.2| , six cases need 
to be distinguished (see Figure ^: 

Case (a), < d < 3/2: As 7 increases from 0, Pt increases smoothly from 1, with an 
initial correction of order 7^, while JF decreases smoothly from 0, with an initial correction of 
order 7^. It follows that ^ is a smooth and decreasing function of Pt in the regime /3r > 1; 
from ( [1.14| ) and ( [4.16[ ) we reobtain the expansion (|0|) as Pt I 1- This of course beats the 
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7 = solution JF = whenever /^^ > 1. On the other hand, for /5t < 1 the only solution 
is given by 7 = and = 0. We therefore have a second-order phase transition at /?t = 1 
(see Figure |^a). [In Section ^.4\ we give explicit formulae for the case d = 1.] 

Case (b), d = 3/2: The parametric plot T{(3t) is qualitatively like that for d < 3/2, but 
since /3r — 1 ~ 7'^ while ~ 7^, the expansion as /?t i 1 has fractional powers ~ — 1)^''^, 
as already seen in ([4.51). There is again a second-order phase transition at /5t = 1 (see 



Figure 0b)- 

Case (c), 3/2<<i<2: A typical parametric plot T{13t) is shown in Figure |^. As 
7 increases from. 0, (3t initially decreases from 1, while T initially increases from 0; the 
dependence T{(3t) is thus initially backward-bending. At 7 = 7cusp, both dPr/dj and 
dj^/d'y pass through zero, so the parametric plot exhibits a cusp. After this, Pt increases 
and JF decreases, with the curve staying below the 7 < •jcusp curve [the latter property follows 
from the monotonicity of the slope ( |3.35| )]; Pt tends to +00 and J-' tends to —00 as 7 t 1/c?- 
Of course, the 7 = solution JF = is also present for all (3t- Since the stable phase is the 
one with the minimal JF, the point where the "cusped branch" crosses the JF = solution 
corresponds to a first-order phase transition with a (i-dependent critical inverse temperature 
Pt,c < 1- This value can be located numerically by solving JF(7) = for 7 using (3.30b), 
and then substituting into (|3.31|) to find Pt- For d slightly greater than 3/2 it is given by 



(E3)/(l 

Cases (d— f), d > 2: These cases are qualitatively similar to case (c). The only difference 
is that for d > 2, reaches a finite value at 7 = l/c?; beyond this point we must consider 



the condensate phase parametrized by / > at 7 = l/c? and g = g^. From ( p.37|) and 
(3.36b) this branch is given by 

T = _/^ + i^gZ±^ (4.18) 

4 2 ^ ' 

(shown as dashed lines in Figures 0d-f). The condensate phase may begin after the cusp 
{d < d ^ 2.55391), at the cusp {d = d) or before the cusp {d > d); this distinction is, however, 
of little interest. What does matter is where the cusped branch crosses the JF = solution: 
for d < d^, ^ 2.38403 this crossing occurs within the "ordinary branch" < 7 < l/rf (Figure 
01), while for d > d^ it occurs within the "condensate branch" f > f^, at •y = 1/d (Figure 0f); 
for d = d^ it occurs precisely at the joining point f = f* (Figure ^). For d < d^,, this crossing 
corresponds to the first-order white-noise-to-isotropic transition; we solve numerically for 
the corresponding (3t,c in the same manner as in case (c). For these dimensions there is also a 
subsequent second-order isotropic-to-nematic transition where the ordinary and condensate 
branches join: this occurs at Pt,c' = fl/^d{f\ — 1)]. For d > d^, the crossing corresponds to 
the first-order white-noise-to-nematic transition; we solve numerically for /9r,c" by solving 
JF(/) = for / using (|4.18| ) [taking only the solution with / > /*] and then substituting into 
( CT V We have, for example, 3t,c"(3) ~ 0.734865, ^t,c"(4) ~ 0.572931, ^t,c"(5) « 0.466744 
and /3t,c"(6) ~ 0.393066. As c/ ^ cx), we have /?r,c" = /^/[2(/oo - 1)]^^"^ + Oid-"^), where 
foo ~ 3.512863 is the unique solution of log / = (/ — l)/2 on / > 1. 
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4.4 One-Dimensional Case 



For the special case ci = 1, we have available the exact formulae ([A.11| )/( |A.12| ) for /(7) 
and g{-f): 

/(7) = il-lY'^' (4.19) 
■1 + (1-72)1/2 

^ ' ^ ' (4.20) 



(4.21) 



^(7) = log 

It is convenient to introduce a new variable s G [0, 1) via 

2s 



7 



/ 



1 + S2 
l + g2 
1 -S2 



(4.22) 



[From ( p?9|) we see that s = e where m is the isovector mass gap.] Using the direct 
approach (|3.12|) , we find 



Hs) = Ms) 



+ 



1 - /3t(1 



,2\2 



-logd-S^ 



The derivative of ( [4. 231 ) with respect to is 

dJ^Q 1 + s' 



d{s' 



2(1 



Pi 



-l + s' 



(4.23) 



(4.24) 



It follows then that for /^^ < 1, JFq is a strictly increasing function of on [0, 1), so the 
absolute minimum lies at = (7 = 0) and has JFq = 0. For /3t > 1, ^0 has a unique local 
minimum at = 1 — 1//?t, which is of course its absolute minimum. We therefore have 



F = inf J^o(s) 

0<s<l 



for < /3t < 1 

-i(^T-l-log^T) for^T>l 



(4.25) 



exhibiting a second-order phase transition at /?t = 1- 

Finally, let us compute the correlation functions. Recalling that in = 1 the free propa- 
gator is a pure exponential 



^-m\x-y\ _ g\x-y\ 



we find from ( p.20|) and ( |2.23| ) the isovector and isotensor correlation functions 



Gvix,y) 



M-y\ 



-.'^\^~y\ 



(4.26) 

(4.27) 
(4.28) 



Of course, the result for Gy reflects the (spurious) spontaneous breaking of Z2 gauge invari- 
ance at (3t > 1, when N = 00. 

In Section ^ we will confirm these results by computing the N ^ 00 limit of the exact 
finite- iV solution. 
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5 Behavior of Solutions for j3y ^ ^ 
5.1 Overview 

We begin by summarizing our conclusions concerning the phase diagram of the N = oo 
mixed isovector /isotensor model. The complete phase diagram is a two-dimensional surface 
in the space^f parameters Pv, Pt, d. In Figures pj-pT^ we show the qualitative phase diagrams 
in the (/^v-, /5T)-plane for different ranges of spatial dimension d: first-order transitions are 
shown in solid lines, and second-order transitions in dashed lines. 

There is always a line of first-order transitions at jSy = 0, /5t > f3T,c{d)- The remaining 
features of the phase diagram depend on the spatial dimension: 

(a) < d < 3/2: There are no phase transitions other than the second-order transition at 
0v, Pt) = (0, 1) and the first-order line at = 0, /5t > 1 (Figure g). 

(b) 3/2 < d < 2: Three curves of first-order transitions meet at the point {i3v,i3t) = 

(0,/9t,c): the first-order line at /3v = 0, /9r > Pt,c', and a pair of symmetrical curves 

Pt,c{,Pv)- Second-order phase transitions occur at the endpoints (/?y, = (±/5\/,e5 /^T.e) 
of the latter curves (Figure 

(c) 2 < d < d^ ^ 2.38403: In addition to the transitions present for 3/2 < d <2, there is 
an additional pair of second-order transition lines located at 

(Figure The intersection of these lines with the Pv = ^ axis at Pt,c' = fll[^d{f^ — 

1)] corresponds to the isotropic-to-nematic transition of the RP^~^ model, while their 
intersections with the /?t = axis at /?y = -^f^/2d correspond to the paramagnetic-to- 
ferromagnetic {Pv > 0) and paramagnetic-to-antiferromagnetic {Pv < 0) transitions of 
the A^-vector model. 

(d) d = d^: The two transitions of the RP^~^ model (at Pt,c and Pt,c') coincide. For 
Pv 7^ 0, the second-order transition lines (|5.1|) remain above the first-order curves 
(Figure 0). 



(e) d^, < d < d^^. = 3: Second-order transitions occur on that part of (|5.1|) where they lie 
above the first-order curves (Figure 0). 

(f) d>?>: The second-order transition lines (|0| ) join onto the endpoints {^Pv^e) PT,e) of 
the first-order transition curves (Figure |T3|). 



Equation (10) of contains an error: where Oku writes K (= our (3t) he should have written K / K'^ 
(= our Pt/ Py); with this change, the equation agrees with our This error can be traced back to a 

notational confusion in an earher paper of the same group | |39| , p. 125], where the authors say "here we have 
scaled a as a ^ iiTa" but without explicitly defining a new symbol; they then presumably confused the scaled 
and unsealed meanings of a. We have been unable to understand where equation (11) of EQ] conies from. 
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The dependence of /?T,e and (5v,e on d is shown in Figures |l^ and |T5[ Together with Figure |T|, 
which illustrates the dependence of (3t,c on d at /Jy = 0, these plots describe the principal 
features of the phase diagram. 

Let us now show how these conclusions can be derived using either the direct or the 
parametric approach. 



5.2 Direct Approach 

We now have to look at the entire free energy ^(7) = ^0(7) + l^v^ii.!) + const. Since 
< foi' all 7, it follows immediately that the absolute minimum of ^^(7) lies in the 
region sgn7 = sgn/5y. In what follows, we restrict attention without loss of generality to 
/3y > and 7 > 0^ 

The plots for /3y = (Figures |^ and ^ get smoothly perturbed as we turn on Pv > 0. 
In particular, for d < 3/2 we observe empirically (though we have been unable to prove it 
analytically) that ^^(7) has a unique local minimum in the region 7 > 0; it is of course the 
absolute minimum. Moreover, it varies smoothly with /?t > and /5y > 0. Therefore, there 
are no phase transitions at /3t 7^ 0. ^ 

For d > 3/2 the story is more complicated. The local minima lying for /5v = at 7 = 
and at 7 = 70 (/?r) both move smoothly as a small /3y > is turned on. Let us call the 
resulting minima 'jaiPv^PT) and 'ybiPv^Pr), respectively. Three nontrivial things can then 
happen: 

(a) A first-order phase transition can occur when the two minima become degenerate 
[^(7a) = ^ilb) with 7a 7^ 7^,]. 

(b) A second-order phase transition can occur when the two minima collide (7a = 7;,) . 

(c) A second-order phase transition to a phase with long-range order can occur when 7^ 
reaches 1/d, provided that it is the absolute minimum [76 = 1/c/ with J-'ijb) < ^{la)]- 

These correspond, respectively, to the first-order curve Pt,c{(^v), the critical endpoints {±Pv,e, /?r,e), 
and the second-order lines I3t,c'{Pv) given by (p.!]). 

One could now construct plots of versus 7 (or versus / for the condensate phase) for 
various values of jSy and Pt- For example, one could look at "slices" at fixed jSy and show 
how the shape of ^^(7) curves is deformed as Pt increases (analogous to Figures || and |^ for 
/3v = 0). One would then see the behavior (a), (b) and (c) corresponding to the first-order 
transition curve, the critical endpoint and the second-order transition line. Rather than 
pursuing this approach in detail, however, we turn to the parametric approach, in which the 
same features can be observed more efficiently because one studies all Pt in a single plot, 
for each fixed Pv- 



5.3 Parametric Approach 

We know that the absolute minimum of the free energy lies in the region sgn7 = sgn/3y- 
We therefore restrict attention to Pv > and consider only 7 > 0. For each fixed Pv, we 
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plot Pt and JF parametrically as a function of 7, using ( p.31| ) and (3.30b). Since fdi'j) is a 
monotonically increasing function of 7 on [0, 1/d], we can (and do) alternatively use / as an 
independent parameter in place of 7; this has the advantage of allowing us to treat normal 
and condensate phases using the same independent parameter / for both of them. We use 
( p.31| )/(3.30b) for the normal phase !</</*, and (|3.37|) /(3.36b) for the condensate phase 
/>/*• 



Figure |I6| contains a "notations guide" illustrating various features of the parametric 
plots of JF versus Pt, to be discussed^presently. In particular, points c ajid c'^correspond, 
respectively, to the first-order curve Pt,c{Pv) and the second-order lines (5t,c'{13v) given by 
(p.ll). The condensate phase f > is shown (both here and in subsequent figures) as a 
dashed hne. 



In Figures p!7|-p^ we plot Pt versus /, and JF versus Pt, at different values of Pv in various 



ranges of dimension d. Two qualitatively different cases need to be distinguished: 

Case (a), < d < 3/2: Here Pt is a strictly increasing function of / on [1, 00) for any 
value of Pv', it tends to +00 as / — > +00 (7 t l/d). Then JF is a smooth function of Pt for 
any Pv 7^ 0. In particular, no phase transition occurs for nonzero Pv- The complete phase 
diagram is schematically shown in Figure (We refrain from showing here the parametric 
plots, as they are so boring.) 

Cases (b— f), d > 3/2: For d > 3/2, the generic behavior Prif) is the following: For 

< Pv < Pv,e there exist /i, /2 G [1, 00) such that Pt is a strictly increasing function of / for 

1 < / < /15 a strictly decreasing function for /i < / < /2, and a strictly increasing function 
for /2 < / < c>o. As Pv increases, the maximum at /i and the minimum at /2 come closer 
together, and they finally collide at / = /e (or equivalently 7 = 7e) when Pv = Pv,e- This 
point is given by the solution of the simultaneous equations dPT/d'j = and d'^PT/d'-y'^ = 
(see Figure ^ for a plot of 7e versus d). Beyond this point, Pt is a strictly increasing function 
of / on the entire domain 1 < / < 00. 

This behavior of PtU) is reflected in the parametric plot J^(Pt), where the values /i and 
/2 correspond to the cusp singularities "1" and "2" of the triangle-shaped structure shown 
schematically in Figure |16[ This is exemplified in Figure 0, which shows the d = 2 case. 
The "triangle" becomes smaller as Pv increases (Figure p!7|a,b,c) and disappears as Pv T Pv,e 
(Figure IT] 



For 3/2<(i<2, a first-order phase transition occurs at the point of self-intersection 
(point c in Figure |l^) for all < Pv < Pv,e- The coordinates Pt,c{Pv) of this transition are 
found by solving the system PT{'ji,d,Pv) = PT{'y2,d,Pv) and J^{'ji,d, Pv) = J^{'j2, d, Pv)- 
They form two symmetric slightly convex branches on the phase diagram, as shown in Figure 
H for the case d = 2.0 The triple point /5t,c(0) is a shghtly decreasing function of d, as shown 
in Figure 0. 

For (i > 2, we need to consider also the branch coming from the condensate phase / > /* 
with 7 = l/(i (shown by a dashed line in Figures T^-P^l). The beginning of this branch 



Magnoli and Ravanini |g, Figure 2] correctly found the endpoint (/3y,e, /3T,e) ~ (0.0134,0.905) with 
7e ~ 0.411, but incorrectly drew the transition curve as concave and as smooth at /3y = 0. They also failed 
to notice the first-order transition curve at f3v = 0, (3t > Pt,c- 
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at point c' (Figure |T6D corresponds to / = A key role is played by the position of c' 
relative to the points c, 2 and 1 (and therefore to the internal structure of the "triangle"); 
this depends on d and on jSy- We have to distinguish several cases: 

Case (c), 2 < d < ~ 2.38403: The evolution of /?t(/) follows the generic pattern 
described earlier for > 3/2. For small values of Pv (Figure 0a,b,c,d) there are two phase 
transitions: a first-order transition at point c of the "triangle" , and a normal-to-condensate 
second-order transition at point c' , which lies outside the "triangle" . These correspond on 



the phase diagram (Figure |T0D to the first-order curve (3t,c{(^v) and the second-order lines 



Pt,c'{i3v), respectively. As we increase the "triangle" shrinks and finally disappears at 



the point {(3v,e, (^T,e) , shown as a dot in Figure [T^. This point corresponds to the endpoint 



of the solid curve in Figure 0, and is a second-order phase transition. Its evolution as a 



function of d is is shown in Figures 14 and 15. For f3v > Pv,e, only the transition at point c' 



remains (dashed line in Figure ITOl) 



Case (d), d = ^ 2.38403: As d increases, the first-order transition curve i3t,c{Pv) 
and the second-order transition lines PT,c'(yPv) come closer together; at d = d^^ they collide 
at Pv = 0, Pt,c = Pt,c' = Pt,* ~ 0.873264 (Figures |1] and |llj). In terms of the parametric 
plot, points c and c' coUide when Py = (Figure |^). For all Pv > 0, point c' continues to 
lie outside the "triangle" , so that the rest of the scenario has the same features as the case 
2 < d < d^, considered above. 

Case (ei), d^ < d < d 2.55391: At small Pv, the point c' lies inside the "triangle", 
between points 2 and c; the first-order phase transition at point c thus corresponds to a 
passage directly into the condensate phase (Figure p!9|a,b). As Pv increases, the condensate 
branch gradually gets "expelled" from the "triangle" ; we define Pv,cc' to be the value of Pv 
at which points c and c' coincide (Figure p!9|c). In Figure ^we plot Pv,cc' as a function of d. 
For Pv > Pv,cc', the point c' lies outside the "triangle" (Figure |19|d), so that the first-order 
transition at point c is a transition within the normal phase (just as in the case 2 < d < d*); 
there is a subsequent second-order transition at point c' to the condensate phase. Finally, at 
Pv = Pv,e, the "triangle" disappears (Figure |I^e); this occurs within the normal phase. For 
Pv > Pv,e the only phase transition is the second-order transition to the condensate phase 



at point c' . The complete phase diagram is shown in Figure 112 



Case (e2), d = d ~ 2.55391: The dimension d is, by definition, the one in which point 
c' coincides with point 2 of the "triangle" when Pv = 0. For any Pv > 0, the description is 
identical to the one given above for the case d^, < d < d. 

Case (es), d <C d <C d** = 3: At Pv = 0, the start of the condensate branch (point c') 
lies between points between 2 and 1 of the "triangle" (Figure pO|a,b). As Pv increases, point 
c' is gradually expelled from the "triangle" : first it reaches point 2 when Pv equals 

Pv,2 = , (5.2) 
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where ( p.39|) vanishes at / = /*; then it reaches point c at Pv = Pv,cc'- Finally, the "triangle" 
disappears at Pv = Pv,e] it does so in the normal phase (since point c' has already been 
expelled from the "triangle"), so that 7e < 1/d (see Figure This evolution is shown in 
detail in Figure pl|a-c for the case d = 2.75. For (3v > Pv,e, the only remaining transition is 
the normal-to-condensate transition at point c' (Figure |20|f). The phase diagram is identical 
to that for d^ < d < d. In Figure ^ we plot l3v,2 < Pv,cc' < Pv,e as a function of d. 



Case (f), d > d** = 3: The three points Pv,2, Pv,cc' and l3v,e all merge at d = d^^ = 3 
[we demonstrate this analytically below). For d > 3, point c' stays inside the "triangle" until 

5). 



the "triangle" disappears at Pv = Pv,e (see Figure ^ for d = 3 and Figure ^ for d 
The phase diagram for > 3 is shown in Figure 0. The dimension d^^, = 3 is equivalently 
characterized by noting that 7e reaches l/rf as | 3 (see Figure |2^) .p^ 

A curious thing happens in dimensions d > 4, where = f'^{l/d) is finite. Then d^T/dj 
for the normal phase, given by (|3.33|), can vanish at j = 1/d: this happens when /3y equals 



k 

2d 



2 - 



(5.3) 



which lies strictly below Pv,2 (since /* > 1 and fl < oo) and corresponds to point c' coinciding 
with point 1 of the "triangle". Indeed, this coincidence occurs for all Pv in the range 
Pv,i < /3v < Pv,2 = Pv,e] it ends only when the "triangle" disappears. See Figures pB|c,d,e 
for an illustration of this phenomenon in c? = 5. 

Let us now prove that d^^ = 3. Using the asymptotic expansion (|A.15|) for fdi'j) as 
7 t l/d, we obtain 

dPj' , ,^ d_i / , ~ N o_d 



where e 



df 

-yd > and 



p,{d) e^-' + p,{d, Pv) e'-^ + 0(e) 



pi{d) 
P2{d,Pv) 



d_Y^ r(l-rf/2) 
2^) d{U - 1) 

27r Uh-dPv) 

d) (/,-i)r(2-rf/2) 



(5.4) 

(5.5) 
(5.6) 



Alternatively, we can write (|5.4| ) as a function of — / as / j 
dPx f*-f ( -27r 



df 



dif. 



1) 



+ P2{d,Pv) 



dT{l-d/2) 



4-d 
d-2 



{h-f)—- + O (/.-/)— . (5.7) 



Let us now fix Pv to the value Pv,2 = /*(2 — f*)/{2d) where point c' coincides with point 
2 of the "triangle" . Kt Pv = Pv,2 the coefficient p2 becomes 



P2{d,Pv,2) 



2tx 



d/2 



(5.^ 



The existence of the special dimension d**, 
by Caracciolo and PeUssetto |l2[. 



2(/,-i)r(2-rf/2) • 

as well as the fact that c?,, = 3, were found independently 
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As e I (i.e. / | /^,), the first term in (|5^ )/ ([5^ ) [which is negative for 2 < (i < 4] dominates 
for d < 3, while the second term (which is positive for < li < 4) dominates for d > 3; 
at d = 3 both terms are of the order e^/^, with net coefficient pi(3) + P2(3, /?y,2(3)) > 0. 
Therefore, for any 2 < c? < 3, we have djST/df < for sufficiently small e > (i.e. for / 
slightly less than /^,). The condition d^r/df < corresponds to the existence of the back- 
bending branch of the "triangle" (between points 1 and 2). It means that the point f]v,e 
where the "triangle" disappears is strictly larger than /3v,2- By contrast, for d > 3, we have 
dPr/df > for sufficiently small e > 0. This means that there is no "triangle" structure for 
any /5y > l3v,2] rather, the "triangle" disappears exactly at /3y = Pv,2- That is, the points 
Pv,2, Pv,cc' and Pv,e all coincide for d> d^^ = 3. 

Finally, as c? — > oo, we have the following asymptotic expressions for the endpoints /3r,e 
and [3v,e'- 



Pv,e 



1 

2d ^ 


1 


7 

8d^ ^ 


15 


153 


2^ ^ 


8# ^ 


32c/5 


1 


1 


3 


33 


3 


2d ~ 


8rf3 " 


8# ~ 


32rf5 


~ d^ 



+ 



451 
32rf6 



+ 



(5.9) 
(5.10) 



5.4 One-Dimensional Case 



The free energy in the direct approach is given by ( |3.12| ). We use the explicit expressions 
( |A.11| ) / ( |A.12|) to eliminate 7, /(7) and ^'(7) in favor of the variable s defined in ( [4.21|) / ( [4.22| ): 



-2s' 



+ 



The derivative of ( p.llD with respect to s is 

(l + s^ 



log(l - s') 



I3v 2s , PI 



+ W- 

/3t 



(5.11) 



dT 

ds 



s + (3t\s + (^](s'-1] 



/3t(1-s2)3 

and therefore the extremum satisfies the cubic equation 



s^ + vs"^ - (1 - r)s - V = 



(5.12) 



(5.13) 



where v = Pv/Pt and r = 1/Pt- For Py > this cubic is strictly convex on s > 0, 
negative at s = and positive at s = 1; it therefore has exactly one root > 0, which 
is simple and lies in < < 1; moreover, s^, is a real-analytic function of Pv,Pt > 0. 
For Py < 0, the relevant solution is obtained by noting that (|5.13|) is invariant under the 



combined substitution P\ 



-Pvi s — i> —s. There is thus no phase transition when Pv 7^ 0. 



Note also that when < /?t < 1, the root s = when f = is simple, so that is a 
real- analytic function of Pv also at Pv = 0. So there is no phase transition when Pv = 
and < Pt < 1- (The same obviously holds also for Pt = and in fact also for Pt < 0.) 
For small Pv > we have: 
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1) For /3t < 1: 



/3- 



V 



Pi 



Pi 



2(1 -/?t) 



(5.14) 
(5.15) 



2) For Pt = 1: 



FiPv,PT) 



\Pv 
3 



i/3_^+ML + o(i/3.n 



|4/3 



+ 



3) For /?T > 1: 



1/2 



i^(/3y,/5T) 



K^T-1 



/9v 



^ 2/3t(/3t - 1) 
log^T) - \Pv\ 



1/2 



0(/5^ 



(5.16) 
(5.17) 

(5.18) 
(5.19) 



We conclude that for /3t > 1 there is a first-order phase transition as Pv passes through 
0: the free energy has an absolute-value cusp, and the isovector energy Ey has a jump 
discontinuity. In particular, at /Jy = the Z2 gauge symmetry is spontaneously broken. Of 
course, this is an artifact of the N —>■ 00 limit, and does not occur for any finite A^. 

Remark. To obtain the free energy in the parametric approach, we plug this into 
either ( p.30a| ) or (3.30b), using the explicit expressions ( |A.11| )/ ([A.12|) to eliminate both 7 
and g in favor of / and hence s: 



-2s2 

1 - S2 



-log(l-s2)+/3Ts' 



loefl 



Pvs 



Finally, by the same reasoning as in Section 4.4 the correlation functions are 



Gv{x,y) 



,2|x-yl 



(5.20a) 
(5.20b) 

(5.21) 
(5.22) 



6 One-Dimensional Case: Pure RP^-\ CP^'^ and QP^'^ 
Models 



In this section we compute the exact finite-iV solution for the one- dimensional RP 



N-l 



CP^ ^ and QP^^ ^ models and discuss the behavior as N 00. We reproduce the N = 00 



jTV-l 
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phase transition at /5t = 1 found in Section |4.4| , and we explain why it is an artifact of the 
N ^ oo hmit. 



6.1 General A^^ 

Note first that the one-hnk integral 







1^ exp 





dVL{(T'2 



(6.i: 



is independent of the unit vector cri [since by G-invariance cri can be rotated to (1, 0, ... , 0)]. 
It follows that for a one-dimensional open chain — or more generally, for any lattice which 
has no closed cycles (i.e. is a forest) — the partition function factorizes as 



(6.2) 



where ^ is the number of links in the lattice. In particular, the infinite-volume free energy 
in one dimension (which is independent of boundary conditions and can thus be computed 
using an open chain) is 



lim 

L— >oo 



--logZ^(/?)^-i 



(6.3) 



So it suffices to study the one-link integral Zisi{[3). 
Let us study the more general integral 



Zn{M) 



dfi(x) 



(6.4) 



where M is an n x n symmetric real matrix and dVt is normalized uniform measure on the unit 
sphere in W^. It is easy to see that Z^iP) equals ZkN^M) where k = 1, 2, 4 for K = M, C, Q 
provided we take M to have k eigenvalues equal to /5/2 and k{N — 1) eigenvalues equal to 0. 



Diagonalizing M and then making the change of variables Xj = x~ in (|6.4|) , we obtain 



exp 



.i=i 



(6.5) 



where /ii, . . . , fin are the eigenvalues of M; here c„ = r(n/2)/27r"'/^ since Zn{0) = 1. Spe- 
cializing now to the case /ii = . . . = /Xfc = /3/2 and fik+i = ... = //„ = 0, we change variables 
to t = Ai + . . . + Afc and obtain 



r(|)r(^) 



(1 -t)- 



dt 



iFi 



k n P 
2' 2' 2 



(6.6a) 
(6.6b) 
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where n = kN and 
iFi(a;c;z) = 



a aia + 1) 

I + -z + —. -f — + 

c c(c+ 1) 2! 



(6.7a) 



Tic) 



1 



Tare 



t) 



c—a—l 



dt 



[for Re c> Re a > 0] 



(6.7b) 



is the confluent hypergeometric function 30, 61|. Thus, 

■1. K. P\ for i?P^-i 
Z7v(/3) = { ^FAl■,N■J^) for CP^-^ 



F (I- K- 

1-^112' 2 ' 2' 

iFi(2;2iV;f) 



(6.^ 



for QP 



N-l 



We remark that the CP 



case can be written in terms of elementary functions: 

N-2 



iFi l;iV; 



(N-iy. 

(/5/2) 



N-l 



E 

i=o 



(/3/2V 



(6.9) 



(A slightly more complicated derivation of the CP^ ^ result has been given by DiVecchia et 
al. 1^ Appendix B.2]. The RP^-^ result is also well known [||, |3|, 3.) 

Finally, let us discuss the correlation functions, restricting attention for simplicity to the 
jlpN-i ixiodel. For the isovector correlation function we have 



Gv{x,y) = (cr^ ■ (Ty) 



5. 



xy 1 



(6.10) 



it vanishes for x 7^ y as a consequence of Z2 gauge invariance. The isotensor correlation 
function can be computed using standard results on one-dimensional 0(A^)-invariant a- 
modelsQ the result is 



GT{x,y) 



{T,-Ty) = ((cr. 



N 



\x-y\ 

N,2 



where 



Vn,2 



P 



l-ri(,2i Y + ^1 9 



' 2' 



1 

A^- 1 
N 



i-' V 2 ' 2 ' 2 ^ 



l-^ll2' Y + -"-i 9 J 



2' 



F (I- S.) 

I-' 1 V 2 ' 2 ' 2' 



- 1 



I-' 1 V 2 ' 2 ' 2' 



+ 



N 



+ 1 



(6.11) 

(6.12a) 
(6.12b) 
(6.12c) 



(the equalities being proven by standard recurrence formulae for iFi] 
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See e.g. g equations (2.33), (3.16) and (3.19)]. 
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6.2 ^ oo Limit 

The large- limit is obtained by taking N oo with j3 = (3/kN fixed.Q For the free 
energy, this limit can be read off by applying Laplace's method to the integral ( |6.6a|) : 



1 

N^'co kN ' ' o<7<l 



F{I3) = - lim —logZN{kN(3) = - sup \[(3t + log{l - 1)] (6.13a) 



for /? < 1 

-i(/3- 1 -log/5) for^>l 



(6.13b) 



in agreement with ( [4.25|) . Thus, in the N = oo limit, the free energy exhibits a second-order 
phase transition at /? = 1: the specific heat Ch ~ d'^F/dp^ has a jump discontinuity. 

On the other hand, we know that such a phase transition is impossible for any finite N. 



Indeed, from the integral (|6.6a]) , it is clear that Zn{(3) is an entire analytic function of /3; 



and since Zm{[3) > for /3 real, it follows that log Zn[I3) is a real-analytic function of 13 (see 
Figure [26|a). How, then, can the nonanalyticity at /? = 1 arise in the limit N oo? The 
answer is that a limit of real-analytic functions need not be real- analytic: singularities can 
come in from the complex plane and pinch the real axis. Since Zj^{(3) is analytic, the only 
singularities of \ogZj^{[3) are at the (complex) zeros of Zjq{(3). So there must be zeros of 
Zn{i3) very close to /3 = 1: 

Theorem 6.1 Fix a > 0. Then there exist, for all sufficiently large c> 0, complex numbers 
z* satisfying i-Fi(a; c; z*) = and lim {zl/c) = 1. 

c—*oo 

In Appendix |B| we give two proofs of this theorem: a "soft" proof based on a direct tran- 



scription of the intuition just sketched (Appendix |B.2| ), and a "hard" proof that provides 

an explicit asymptotic expansion for z* in powers of c~^/^ (Appendix |B.3| ). Surprisingly, we 

have been unable to find this result anywhere in the extensive mathematical literature on 

confluent hypergeometric functions. 

For the isotensor correlation function (|6.11|) / (|6.12|) , we can apply Theorem p.l| to (|6.12c| ) 

to obtain ^ 

/ txtTTn I for /3 < 1 

hm vnA^P) = < ~ ~ - 6.14 

[1-1//3 ioT/3>l 

in agreement with (^4.28|) . The approach to the N ^ oo limit is shown in Figure p^b. 



7 One-Dimensional Case: Mixed Isovector/Isotensor 
Models 

Let us now generalize these one- dimensional results to the mixed isovector /isotensor 
model (0). 

■^^ Some authors define (3 — P/N . But the present definition is more suited to treating the RP^^^, CP^~^ 
and QP^~^ cases in a unified manner. 
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7.1 General N 

Consider the mixed isovector/isotensor Hamiltonian 



a I ^2 



{7.1] 



with (T^ G ^, on a one-dimensional lattice of L sites with free boundary conditions. We 
obtain immediately the partition function 



where £ = L — 1 is the number of links and 



exp 



Pi 



Pv(Ti ■ (72 + ^(CTI ■ 0-2)^ 



is the one-link integral. The free energy is 



Fn{Pv,Pt) = lim 



L—i-oo 



-^\0gZN{Pv,PT)''-' 



- log Zn{Pv,Pt) 



(7.2) 



(7.3) 



(7.4) 



To evaluate the one-link integral ( [7. 3D , let us fix cri = (1,0,..., 0) and write (T2 
[si, . . . , sn), so that 



Zn{Pv, Pi 



T{N/2) 
27r^/2 



r(f 



r(i)r(^) 



j S{sl + ■ ■ ■ + 4 - 1) e/^^^i+(/^^/2)^? ds,--- ds 

(l_s2)(7V-3)/2g^p 



In the two limiting cases /?t = and Pv = 0, we have 



z^{Pv,Q) = r(iv/2) ( — 



(N~2)/2 



Pvs + 



J{N-2)/2iPv 



ds . 



Zn{0,Pt) = iFi 



l_N_pT 
2' T' ~2 



(7.5a) 
(7.5b) 

(7.6) 
(7.7) 



where Ii, is a modified Bessel function. Expanding (7.5b) in powers of either Pt or Pv, we 
get the series representations 



and 



Zn{Pv, Pi 



Zn{Pv, Pi 



^{Pt/2)^ (1/2), ^(f^.i.i.N PI 



k=0 



k\ iN/2)k 



fc=0 
00 

E 

fc=0 



(2A;)! (iV/2)fc 
(/5y/2)2'= 1 



k\ (iV/2), 



2 ' 2 



(7.8) 

(7.9a) 
(7.9b) 
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where (a) 



r{a + n)/r{a) is Pochhammer's symbol. Or we can have the double series 

(1/2),+, i(3Tm/3'y/Ay 



Zn{(3v-, Pi 



E 

kJ=0 



{N/2)k+i{l/2)i 



kin 



(7.10) 



which is a special case of the Kampe de Feriet double hypergeometric series p. 150, 
equation 29] ^ Section 1.5] ||, pp. 26-27]. 

It is worth remarking that Z]\r{(3v, Pt) has a kind of "Lee- Yang property": for all > 1 
(not necessarily integer) and all Pt > 0, the zeros of Z^i^Pvi Pt) in the complex /3y-plane 
all lie on the imaginary axis. Indeed, for Pt = Q this is a well-known property of the Bessel 
functions ( |7.6| )p^ while for /5r > it follows from the formula 



Zn{Pv,P 



T 



exp 



V 2 dPl. 



2n{Pv,0) 



(7.11) 



combined with the fact that the operator exp[(/3T/2) (d'^/dPy)] preserves the Lee- Yang prop- 
erty This Lee- Yang property fails for sufficiently negative Pt, except when = 1 



Let us now compute the isovector and isotensor correlation functions. Standard results 
on one- dimensional 0(A^)-invariant cr-models yield 



Gv{x,y) 



{(Tx- O- y) 



,\^-y\ 



((cr^ ■ (Tyf) - — 



,\^-y\ 



(7.12a) 
(7.12b) 



where 



Vna 



1 

/ 

-1 



Js(l_s2)(iV-3)/2 g/3v..+(/3W2).2 



J (1 _ g2YN-3)/2 g/3vs+(/3T/2)s2 
-1 



(7.13a) 



Vn,2 



J 1^ 



/N 



j2^(Af-3)/2 g/3vs+(/3T/2)s2 



J (1 _ _52)(Af-3)/2 g/3t,s+(/3T/2)s2 ^g 



(7.13b) 



For the A^-vector model [Pt = 0) we recover the well-known result 



VN,k 



t{N/2)+k-l{Pv 
I{N/2)-l{Pv) 



while for the pure RP model {Pv = 0) we recover (|6.12| ). 



(7.14) 



'^^ The zeros of lu are pure imaginary for all v > —1 section 15.27], so the right-hand side of (7.6) has 
the Lee- Yang property for all > 0. But the integral (7.5b) is convergent only for > 1 (or, by a limiting 
process, also for iV — 1). 
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7.2 ^ oo Limit in Mixed Models with Various Scalings 

Let us now consider the limit ^ oo with Pv and Pt scaled as 

py = N'py (7.15a) 
P^ = NPt (7.15b) 

where < e < 1. We apply Laplace's method to the integral (7.5b), which can be written 
in the form 



-2' 



Zn{Pv,Pt) = ^n. 'Li, I {l-sy'^'e-'^^^^Us (7.16) 



where 



^(.) = log(l - s') - ^s' - N^-'Pys . (7.17) 



If < e < 1, the behavior at = oo is exactly as in the pure RP'^ ^ model {Pv = 0): 
for /9t < 1, ^{s) has an absolute minimum at s = 0; while for Pt > 1, ^{s) has a local 
maximum at s = and absolute minima at s = ±so, where 

1 \V2 



So = 1 - ^ . (7.18) 

V Pt/ 

There is therefore a second-order phase transition for Pt = 1, and the value of Pv plays no 
role. 

However, for e = 1 the situation is very different. The stationary points of the function 
jF(s) are the solutions of the cubic equation 

+ vs^ - {1 -r)s -V = , (7.19) 

where v = Pv/Pt and r = 1/Pt- The absolute minimum of JF is given by the unique root Sj, 
of (|7.19|) lying in the interval < (sgn/9y)s < 1.0 This is exactly the same cubic equation 
( p.l3| ) found in Section even though the free energy is different [compare ( |5.11 ) with 
( [7.17] )]. PI We therefore draw the same conclusions: there is no phase transition when Pv 
or < Pt < I', but for Pt > I (resp. Pt = 1) there is a first-order (resp. second-order) phase 
transition as Pv passes through 0. ^ ^ 

The phase transition observed as Pv passes through when /3r > 1 can be understood 
as arising from pure- imaginary zeros of Zjy^Pv, Pt) in the complex /3y-plane that approach 
the origin as ^ oo. Indeed, Laplace's method apphed to (7.5b) immediately yields 

hm ^^(/^^'^/^) = cosHPvSo) (7.20) 
Z^{0,NPt) 
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Using elementary monotonicity arguments one can show that for any (3t and nonzero (3v, equation 



( 7.19 ) always has a unique real root in the interval < {sgn f3v)s < 1. 
27 We remark that (|7l7|) = 2 (5. 20b) - K20^. 
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where _ 

for (3t < 1 

^0 = { , , xi/2 ^ (7.21) 

for pT > 1 



1 



It follows that if /3t > 1; 2Ar(/3y, N^t) has (for all sufficiently large A^) zeros in the complex 
/3y-plane that tend to ±{2k + l)7ri/2so as — oo. (As mentioned previously, we know that 
these zeros lie on the imaginary axis for all A^.) In the variable /3y = (3v/N, these zeros tend 
to the origin at a rate 1/A^, producing a phase transition in the limit. For < /?r < 1, by 
contrast, the zeros in the /Jy-plane scale with A^ and thus stay a finite distance away from 
the origin in the /5y-plane. Indeed, the closest zeros to the origin lie on the imaginary axis 
at the Pv value where two saddle points given by (|7.19|) collide, namely 



~ (r - 4 + Vr^ + 8r) y^r + 2 - Vr^ + 8r 

= ± ^^-''^ 

where r = 1/(3t > 1; this goes monotonically from ±i/2 at /^y = to zero at /3t = 1 
(see Figure |2^).0 The corrections to (|7.22|) are presumably of order N~'^^^, as they are at 
/3t = For /3t = 1, the zeros in the /Sy-plane presumably scale as a fractional power of 
A^; empirically we find that they scale like N^^^, exactly like the leading correction term for 
< /3r < 1 (see Figure pS]). 

Let us conclude by considering the correlation functions (7.12) / (7.13) in the limit N oo 
with PvjPt fixed. In this limit, the integrals in both numerator and denominator of (7.13) 
are dominated by s*, so we obtain 



N 



lim VN^k{N(3) = si (7.23) 



for /c = 1, 2, in complete agreement with ( [5. 211 ) / ( [5.221 ). We thus verify by explicit calculation 



the factorization (|2.23| ) in the N —>■ oo limit. Of course, this solution exhibits the (spurious) 
spontaneous breaking of Z2 gauge invariance at /3y = and /3t > 1, when N = 00. 



8 Conclusions 

What is the phase diagram of the mixed isovector/isotensor model ( |1.2D for large but 
finite A^? We know some things for sure: 

With a bit more work one can presumably show that the free energy F 



■YniLN->oo{'^/N)\ogZN{N Pv tN Pt) has a branch cut on the imaginary axis starting at the points (7.22) 



From this it would follow, by the methods of Appendix |B.2| , that zeros of Zn{NPv , N(3t) accumulate densely 
on this cut as iV — s- cx). 

At /3t = we have ZAr(/3y, 0) '-^ I(n~2)/2[Pv), and the first root of Iv{z) is given hy z = where 
~ + 1.855751/1/3 + 1.03315i/"i/3 - 0.003971/"^ + ... 

ijssl , section 15.8] equation 9.5.14]. We have no proof that this same scaling of the correction terms 
persists for < /3t < 1, but it seems plausible and is consistent with our numerical calculations. 
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(a) In dimension li < 1, no phase transitions of any kind are possible in the finite {/Sv, Pt)- 
plane [||, 



(b) In dimension < 2, no spontaneous breaking of the SO{N) global symmetry is possible 
in the finite /3T)-plane |^ |70|, |7ll. In particular, magnetic and nematic long-range 
order are impossible, and all correlation functions decay at least as an inverse power 
of the distance [|73, • 

is possible in any 



(c) No spontaneous breaking of the Z2 gauge invariance at Pv 
dimension 



In addition, some properties of the phase diagram in dimension d > 2 are known rigorously: 



(d) All correlation functions ((cr^^ 



Vn 



)) in the mixed isovector/isotensor 



model at {(3v,Pt) are bounded above by the corresponding Ising-model correlation 
functions {exj^ey-^ ■ ■ ■ £xn^y„) at /5 = \Pv\- we prove this in Appendix |C| (Theorem |C5 



and Corollary |C.6| ). In particular, for \Pv\ ^ Pcising there is no ferromagnetic or 
antiferromagnetic long-range order. 



(e) For Pv > Nf^:/{2d) there is ferromagnetic long-range order 
Pv < —Nf^/{2d) there is antiferromagnetic long-range order. 



Likewise, for 



(f) For (3t > ^2iN~i)^ M there is nematic long-range order Q pi 



Facts (e) and (f) are proven by the method of infrared bounds |7^, |7^, It is noteworthy 
that (e) is asymptotically sharp as ^ 00 when /3t = 0: it captures the leading large- A^ 
behavior ( ^.1|) of the pure A^-vector model. By contrast, (f) is a very poor bound for large 
A^: it behaves as A^^ while the true critical point is presumably of order A^. 

One of the most striking aspects of the N = 00 solution is that it violates fact (c): 
spontaneous breaking of the Z2 gauge invariance occurs at A^ = 00 whenever the absolute 
minimum of the free energy in the pure RP^~^ model is located at 7 7^ [cf. the discussion 
surrounding ( p.l7| )]. But the distinction at A^ = 00 between the white-noise phase (7 = 0) 
and the isotropic phase (0 < 7 < l/rf) does not survive to finite A^. Indeed, at finite A^ 
the correlation functions never have the white-noise form (except of course at I3t = 0). 
Rather, the white- noise "phase" Q < (3t < Pt,c (or < Pt < Pt,c" in dimensions d > d^) 
simply corresponds to the region where 7 is of order 1/A^ [and hence the mass gap ( |3.9| ) 
is of order logA^] when A^ ^ 00 at fixed (3ti while the isotropic "phase" corresponds to 
the region where 7 and the mass gap are of order 1 when A^ ^ 00 at fixed (3t- There is 
no sharp distinction between these two regions at any finite A^; rather, for each finite A^, 
the white-noise and isotropic "phases" both belong to a single paramagnetic phase in which 
the isotensor mass gap itlt is finite and nonzero. This behavior is seen explicitly in the 
d = 1 solution (|6.1 1|) / (|6. 12|) / ( |6TT4D [see Figure p^]; we conjecture that it constitutes the 
correct interpretation of the "white-noise-to-isotropic (pseudo)transition" in all dimensions. 
Likewise, the spontaneous breaking of the Z2 gauge invariance in the isotropic "phase" at 
N = 00, arising from the absolute- value cusp in the free energy F = Fq — 0{\Pv\) as Pv 
passes through 0, cannot survive to finite A^; the behavior at small Pv must be smooth [^ . 
Rather, we can understand the white- noise "phase" < Pt < Pt,c (or < Pt < Pt,c") as 
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corresponding to the region where this rounding takes place over a (3v interval of order A^, 
while the isotropic "phase" Pt > Pt,c (or jSr > Pt^c") corresponds to the region where this 
rounding takes place over a (3v interval of order o{N). Once again, this behavior is seen 
explicitly in the d = 1 solution, where the rounding in the isotropic "phase" takes place over 
a Pv interval of order 1 [cf. the discussion surrounding (|7.20|) /( [7.21|) ]; we conjecture that this 



constitutes the correct interpretation of the "spontaneous breaking of Z2 gauge invariance" 
in all dimensions. 

Putting facts (a)-(f ) together with our analysis of the large- limit — and its pathologies 



in d = 1 — we are led to conjecture the phase diagram shown in Figure for the finite- 
mixed isovector/isotensor model in dimensions d > 2. When both Pt and \Pv\ are small, 
there is a paramagnetic phase in which both the isovector and isotensor correlation functions 
decay exponentially; in particular, there is no long-range order of any kind. Along the Pt = 
axis (pure A^- vector model), there is a second-order paramagnetic-to-ferromagnetic transi- 
tion at Pv = Pc,N-vector, and a second-order paramagnetic-to-antiferromagnetic transition 
at Pv = —Pc,N -vector', foT large A^ we have Pc,N-vector = N f^./{2d) + 0(1). At the opposite 
extreme Pt = +00, the spins are forced to point up or down along a single axis, so the 
0(A^)-invariant correlations are those of an Ising model with coupling Pv'-, therefore, there is 
a second-order paramagnetic-to-ferromagnetic transition at Pv = Poising and a second-order 
paramagnetic-to-antiferromagnetic transition at Pv = —Pcjsing (independent of A^). It is 
reasonable to conjecture that the paramagnetic-to-ferromagnetic transitions at /3t = and 
Pt = +00 are the endpoints of a continuous curve of paramagnetic-to-ferromagnetic transi- 
tions; moreover, this curve Pv,c{Pt) should be monotone decreasing (since Pv and Pt both 
promote ferromagnetic ordering). In the entire region to the right of this curve, there is 
ferromagnetic long-range order. It is also reasonable to conjecture, though this is much less 
certain, that the transition remains second-order everywhere along this curve. There is, of 
course, a symmetrical paramagnetic-to-antiferromagnetic transition curve in the half-plane 
Pv < 0. 

Finally, along the /3y = axis (pure RP^~^ model), there is a paramagnetic-to-nematic 



transition at Pt = Pc^rpn~i. On the basis of mean-field theory [g^, |5^, Q and Monte 
Carlo simulations (for = 3, A^ = 3, 4, 5) this transition is expected to be first-order 
for all A^ > 2, at least in dimensions d > 3.0 (Perhaps the transition is second-order for 
2 < d < (i*.^ We expect that P^^rpn-i scales for large A^ proportional to A^; we predict 
(but with less confidence) that the limiting proportionality constant is Pt,c' [defined by ( |4.9|) ] 
when 2 < d < d^,, and Pt,c" [defined by ( [4.13|) ] when d> d^. 

What about the first-order transition curve Pt,c{.Pv) found at A^ = cxd in dimensions 
d > 3/2 (Figures P|-[TB|) — does it survive to finite A^? We doubt it: this transition is 
the extension to Pv ^ of the "white-noise-to-isotropic transition" of the RP^~^ model, 
which we have argued is an artifact of the N = 00 limit. Moreover, extensive Monte Carlo 



■^^ By contrast, for N = 2 the RP^ ^ model is equivalent, via the change of variables 9' = 29, to the 
standard XY model, so the transition is second-order. 

■^^ Of course, it is far from clear what noninteger dimensions d should be taken to mean nonperturbatively. 
But one can also study the 2-dimensional model with power-law-decaying long-range interaction K^y 
\x — (cr > small), which should behave qualitatively like the putative short-range model in 

dimension d = 2 + a. 
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simulations of the two-dimensional mixed isovector/isotensor model for = 3, along with 
preliminary simulations of the two-dimensional RP^~^ models for = 4, 6, 8, showed no 



sign of any first-order transition |]16 



One plausible scenario is the following: All transitions (whether first-order or second- 
order) within the "normal phase" I7I < 1/d are artifacts of the N = 00 limit, but transitions 
(whether first-order or second-order) into the "condensate phase" I7I = 1/d do survive to 
finite A^. In this scenario, the transition curve in the limit N 00 would be given by 
the upper envelope of the transitions shown in Figures |Tn|-0: namely, by the line i3t,c'{Pv) 
defined by ( p.!] ) for 2 < d < d^, and by max[/3T,c(/3y), Pt,c'{Pv)] ioi d > d^,. Most likely, this 
curve corresponds to the part TF-Fl of the paramagnetic-to-ferromagnetic transition curve. 
In particular, point TF would lie at a Pv value that is o{N) [e.g. most likely of order 1] and a 
Pt value that is close to N(3t,c' [resp. N(3t,c"] for 2 < d < [resp. d > d^], so that as A^ cxd 
the whole curve N-TF is compressed to the point Pv = 0, Pt = Pt,c' [resp. Pt = Pt,c"]- 
Moreover, the whole nematic region (bounded by AF2-TA-N-TF-F2) is compressed at 
A^ = 00 to the line Pv = 0, Pt > Pt,c' [resp. Pt > Pt,c"]', this compression is what produces 
the spurious first-order transition as this line is crossed. It would be useful to test these 
predictions by Monte Carlo simulation of the 3-dimensional mixed isovector/isotensor model 
for various values of A^. 

In dimensions < 2 we expect that there are no phase transitions of any kind in the 
finite (/^v, /5T)-plane. However, for 1 < < 2 there is an Ising-like phase transition at 
Pt = +00, Pv = Pcising] and it is at present unclear whether "new" continuum limits — 
that is, continuum theories whose isotensor correlation functions differ from those of the 
continuum limit of the pure A^-vector model — can be obtained by approaching this point 
|T31 (see also 



ni 



for the analogous problem in d = 1). 



A The Lattice Propagator and Free Energy 



A.l Formulae 

In this appendix we record some needed formulae concerning the lattice propagator and 
free energy 



9d{l) 



d'^p 
d^p 



-\ -1 



1 - 7^ cos Pi 



i=l 



log 



1 — 7 cos Pi 



i=l 



(A.l) 



(A.2) 



which are well-defined for I7I < 1/d [with /^(il/fi) = 00 for d < 2]. Inserting the integral 
representations 



00 



e-'' dt 



(A.3) 



38 



dt (A.4) 



and performing the momentum-space integrals, we obtain 



oo 



/,(7) = le-'lo{itYdt (A.5) 





oo 





where /q is a modified Bessel function. This defines an analytic continuation to non-integer 
d. 



Note also the relations 

dUl) 



dj 



dgdii) fdii) - 1 



d'j 7 
which can be obtained by writing 



- [{t-l)e-'lo{it)Ut (A.7) 
7 J 

(A.8) 



and integrating by parts. In particular, we have 

oo 

7/d(7) + /d(7) = [ te-' loiit)" dt 





oo oo 



j dh j dt2e-''-'' Io(-f{h + t2: 





oo oo 



> Jdt^J rft2e-*^-*^/o(7ti)'/o(7t2)' 




\2 



= UiY (A.io) 

for all (i > (with strict inequality for d > 0), since logJo is a convex function and satisfies 
log/o(0) = 0.0 It follows that (/^ — l)/7/d is an increasing function of 7. 
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The convexity of log/o follows easily from Holder's inequality applied to the integral representation 

7r 
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The cases d = 1,2 can be expressed in terms of standard functions: 

;i_^2)-i/2 forc; = i 

^K(27) for d 



where is a complete elliptic integral of the first kind. The case d = 3 can also be expressed 
in terms of the elliptic integral K, but the result is a god-awful mess Similarly 
we have0 



logfii^^l ford 



9dh) = { \ ' J (A.12) 

-^4^3(1, |, 1,1; 2, 2, 2; 47^) for d = 2 

In general dimension we have the expansions for small 7: 

d . 3d(2d-l) , 5d(6c/2-M + 4) . 

(7) = 1 + -7 H ^ -1 H ^ -1 

\u 2' ^ 8 ' 16 ' 

35rf(24rf3 _ 72c/2 + 82rf - 33) § 63rf(120rf^ - 600c/3 + 1250^^ - 1225c/ + 456) 

^ 128 ^ ^ 256 ^ 

231c/(720c/^ - 5400c/^ + 17700^^ - 30600^^ + 27041c/ - 9460) 

^ 1024 ^ 

42M(5040ci^ - 52920c/5 + 249900c/^ - 661500c/3 + 1011017c/2 - 826336ci + 274800) 

H 7 + 

2048 ' 



(A 



c/ 2 3c/(2c/-l) 4 5c/(6c/2 - 9c/ + 4) g 
(7) - - ^7 7 7 

35c/(24c/3 - 72c/2 + 82ci - 33) g 63c/(120c/^ - 600ci3 + 1250c/2 - 1225ci + 456) 

1024 ^ 2560 ^ 

77d{720d^ - 5400c/^ + 17700ci3 - 30600c/2 + 27041c/ - 9460) ^2 

4096 ^ 
429c/(5040c/6 - 52920c/5 + 249900ci^ - 661500ci3 + 1011017ci2 - 826336c/ + 274800) 

28672 



7" 



(A 



We remark that for integer d > 1, the Taylor coefficients of fd and —gd are all positive, as can 
be seen by expanding ( [A.ID and ( |A.2| ). But our numerical computations of the expansions 



( |A.13| )/( |A.14| ) up to order 7^^ suggest that this property holds only for integer d (that is. 



for every noninteger d, at least one of the Taylor coefficients is negative). 

(see e. g. [fz^ , Theorem 1.2.3]). Note also that for integer d> 1, the inequality ( |A.10[) can be proven directly 



from ( A.l ), using the Schwarz inequality. 

■^•^ The explicit expression for 32(7) is obtained by taking the limit .92(7) = lime^o [r(e) ^ 7)], 
where R{e,j) = t^~^e~H§ijt) dt = r(e)4F3(i, |, f , ^; 1, 1, 1; 47^) = r(e)F4(f , ^; 1, 1; 7'; 7') = 

T^kLo fc!r(fc+i) 2^i("fe'"fe; 1' l)(7/2)^^ iol. Here Fi{a,b; c, c'; z, ^) is an AppeU's function. Since 

2Fi(-fc,-fc;l,l) = r(2fc + l)/r2(fc + 1), we obtain 52(7) = -EZi ^^'fJ?'''^ i^/^)'' = 

— i J^^ 3^2(1, |; 2, 2; x)dx = — ^4F3(1, |, |, 1; 2, 2, 2; 47^). An alternate expression in terms of 3^2 can be 
deduced from equation 18]; we thank Professor Glasser for drawing our attention to this. 
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As 7 T we have the asymptotic expansions 

' {d/2TTy/^T{l - f )(1 - rf7)('^-2)/2 ^ 0(1^ for < d < 2 

fdil) = i -ilog(l-27) + f log2 + o((l-27)log(l-27)) for = 2 (A.15) 

^ + {d/27TY/^T{l - f )(1 - rf7)('^-2)/2 + 0(1 - rf7) for 2 < d < 4 

The hmiting value = /^(l/d) is infinite for < 2 and finite for d > 2. It has the 

known exact value for d = 3 [79|n 

12 

/*(3) = — (18 + 12V2-10v^-7V6)ir^[(2- v^)(V3- v^)] (A.16a) 

r(^)r(A)r(^)r(ii) (A.ieb) 



327r3 
73- 1 



r^(^)r^(i) (A.16C) 



327r3 

^ 1.51638 6059151978 (A.16d) 

and the following approximate values for higher dimensions: 

/*(4) ^ 1.23946 71218 (A.17) 
/*(5) ^ 1.15630 81248 (A.18) 
/,(6) ^ 1.11696 33732 (A.19) 

As d ^ oo we have the asymptotic expansion 

, , „ 1 3 3 15 355 595 , , , 

The limiting slope fl{d) = lim dfdi'y) / d'j is infinite for < 4 and finite for d > 4. Its 

in/d 

asymptotic expansion as — > oo is 

3 39 285 2325 , 

= i + 5 + 4d5 + SP + 

The limiting second derivative f'J{d) = lim d"^ fj_{'j) / d'j'^ is infinite for c/ < 6 and finite for 

7Tl/d 

(i > 6. Its asymptotic expansion as (i — cx) is 

= ^ + 9 + ^ + ^ + ■ (A-22) 



■^'^ The formula in terms of the eUiptic integral is due to Watson ||8^ ; the formula in terms of four gamma 
functions is due to Glasser and Zucker ||79| ; the reduction to a formula involving only two gamma functions 
is due to Borwein and Zucker js^. See also for recent extensions. (N.B.: The original paper ||7^ has 
a typographical error in the numerical prefactor, as does [^4] , p. 21]. The correct prefactor, given here in 
(A. 16b), has been reported in Issl, p. 126] and ]§^, erratum].) 
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The limiting value g*{d) = gd{l/d) is finite in all dimensions d > 0. It has the exact 
values 



^=,(1) = -log2 -0.693147 (A.23) 
AG 

g^{2) = 21og2 ^ -0.220051 (A.24) 

TT 

where G ^ 0.915966 is Catalan's constant, and the following approximate values for higher 
dimensions: 

^,(3) ^ -0.1183701593 (A.25) 

^^(4) ^ -0.07973 38971 (A.26) 

^,(5) ^ -0.06009 70332 (A.27) 

g^{6) « -0.0482796878 (A.28) 

Finally, as o? — > oo we have the asymptotic expansion 

_ 3 7_ _ 45 _ 269 _ 805 

^ ~ M 16^2 32# 128^4 384^5 512^6 + •• " ^^'^^^ 

A. 2 Numerical Calculation of fdil) and its Derivatives 

In this section we discuss briefly how to make accurate numerical computations of fdi'j) 
and its derivatives. What is slightly tricky is to maintain accuracy when |7| is near 1/d. 

The function fdil) and its first and second derivatives can be written in the following 
convenient way: 

/,(7) = Ao{j,d) (A.30) 
7/^(7) = A,ij,d)-Ao{j,d) (A.31) 
7V;(7) = A,{j,d)-4A,{j,d) + 2Ao{j,d) (A.32) 

where 

oo 

^fc(7,d) = j t^e-'h{-itYdt . (A.33) 



For large values of t, the integrand behaves as t^e.~^'^~'^d)t ^ ^^laX the integral converges slowly 
when 7 is close to 1/d. To improve the convergence, we split the interval of t-integration 
into two subintervals [0, a] and [a, oo), where a is chosen in such a way that for z > a/d the 
function Io{zY is well approximated by its large- 2; asymptotic expansion 



dz 



(27r^)<^/2 



d d{d + 8) d{d'^ + 2Ad + 200) 
Sz ^ 128^2 + 3072^3 



+ 0(1//) 



(A.34) 



Introducing the subtracted quantities 



Fo(z,d) = ^o(^)'^- (A.35) 



42 



dz / J 



we can then rewrite ( |A.33| ) as 

Ao{l,d) = [ e-*/o(7t)'rft + / e-'Fo{jt,d)dt + iin^^^^ r(l - f , a(l - 7^)) 
J J 27r7 2 

a 

(A.38) 

Ai(7,rf) = [ te~^Io{-ftydt + I e'HFi{-ft, d) dt + ~ ^^^^ ^(2 - f , a(l - -fd)) 

a 
a 00 , 

r r (i-'jd)^'^ 

A2{j,d) = e-Vlo{it)Ut + e-H^F2{jt,d)dt + ^ r(3 - f , a(l - 7^^)) 

J J (27r7) 2 

a 



where 

T{a,z) = j e^^e'^dt [AAl] 

z 

is the incomplete gamma-function. In particular, at 7 = l/rf we have 

a 00 4 

fM = fd{l/d) = j e''h{t/dYdt + j e-'Fo{t/d,d)dt + j^(^-^y a'~l (A.42) 

a 

for d > 2. Using the above formulae, one can compute /*(3) ~ 1.5163860591520 and compare 
it with the exact value (|A.16|) . 

Similar techniques can be used for Qdi'j), but they are not really needed unless d is very 
close to zero. 

Remark. An alternative method for computing these and similar quantities to very high 
precision, with rigorous error bounds, is given in ||8^, Appendix B]. See also Appendix 
A and erratum] for a brief summary. 
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B Confluent Hypergeometric Function iFi(a; 6; b^) in the 
Limit 6 ^ 00 



B.l Asymptotic Expansion in Powers of 1/6 at Fixed ^ 

In this appendix we shall develop an asymptotic expansion for the confluent hypergeo- 
metric function iFi{a;b;b^) in the limit b +00, where a is a fixed positive real number 
and ^ = X + iy is a fixed complex number. Our analysis will be based on the integral 
representation 

1 

iFUa-b-z) = ^-r^}^ r [ e'' (1 - tf-^-Ut [for Re&> Rea > 0] (B.l) 

T{a)T{b - a) J 



combined with Laplace's method for contour integrals |8^, sections 4.6-4.10]. We have 



iFi(a; b; b^) = C{a, b) j e^^^ ^a-i (1 _ , (B.2) 



where 

C{a,b) = (B.3) 
i (a)i (0 — a) 

F{t) = et + log(l-t) (B.4) 

The path of integration in ( p.2|) runs initially along the straight line segment from t = to 
t = 1, but we are free to deform the contour of integration in the complex t-plane. Writing 
t = u + iv, we have 

ReF{u,v) = xu — yv + ^\og[{l — u)'^ + v"^] (B.5a) 

lmF{u,v) = xf + — arctan^^^ ^ (B.5b) 

In particular, at the endpoints we have F{0) = and Re-F(l) = —00. There is a unique 
saddle point F'{t^) = at t^, = 1 — 1/^, i.e. 

= l-x/{x^ + y^) (B.6a) 

V. = y/{x^ + y^) (B.6b) 

The value of F at this saddle point is 

F{U) = e-l-loge, (B.7) 

and we have in particular 

Re F{Q = x-\-\ log(a;2 + y'^) . (B.8) 
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For reasons that we will explain shortly, the curve ReF(tj,) = plays a special role in 
characterizing the asymptotic behavior of i-Fi(a; b; 6^) in the limit b oo. Let us therefore 
denote by iS_ (resp. 5+) the part of the curve 



X 



1 - ilog(x2 + = (B.9) 



that lies in the region x < 1 (resp. x > 1), and write 5 = iS_ U 5+ (see Figure To the 

left of i.e. in the region 

A = = X + iy: X < 1 or X - 1 - J log{x^ + y^) < 0} , (B.IO) 

we shall show that i-Fi(a; b; 6^) has a finite limit as 6 ^ oo. To the right of 5+, i.e. in the 
region 

B = = X + iy: X > 1 and x - 1 - | log(x^ + y"^) > 0} , (B.ll) 

we shall show that iFi{a; b; bC,) grows exponentially as 6 ^ oo. More precisely, we shall prove 
the following asymptotic expansions: 

Theorem B.l 

(a) For G A, we have the asymptotic expansion 

a{a + l){a + 2){a + 3) f ^ V ^ 2a(a + l)(a + 2) / ^ V ^ a{a + V , " ^ 



+ 0{b-'')\ (B. 



uniformly for ^ in compact subsets of A, where (1 — ^) is defined as the principal 
branch (i.e. arg(l — ^) = for real < I). 

(b ) For ^ & B, we have the asymptotic expansion 



,F,{a-b-M) = e''(«-i-'°^«)6'^^^^^(e-ir^ x 

r(a) 



126 V ^-1 (^-1)^ y ^ V / ; 

(B.13) 

uniformly for ^ in compact subsets of B, where log^ and {C, — ly^^ are defined as the 
principal branch (i.e. arg^ = arg(^ — 1) = for ^ real > I). 



■^^ This curve was introduced by Szego in connection with the zeros of the partial sums of the Taylor 
series of [that is, Sn(z) — J2k=o 1^^] £^nd of the corresponding remainders e^ — Sn(z). See also |9^, p2[ . 
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Let us first consider region A, which it is convenient to subdivide as follows: 



^1 
As 



{^ = x + iy 
= x + iy 
= x + iy 



X < 1 and X — 1 — I log(x^ + y"^) 



>0} 



a; < 1 and x — 1 — ^ log(x^ + y^) < 0} 
a; > 1 and x — l — \ log(a;^ + y^) < 0} 



(B.14a) 
(B.14b) 
(B.14c) 



Consider first the region Re^ < 1. It is easy to see that ReF(t) < for all real t G (0, 1] 
so the endpoint t = is the dominant point on the line segment of integration in ( [B.2| ) 



By standard theorems on Laplace- type integrals |^ section 4.6], we obtain (|B.12|) provided 
that ^ ^ 1. This proves ( [BTT^ for ^ e ^1 U 

Let us now demonstrate that this expansion is valid in the entire region A. If Re^ > 1, 
the real line segment < t < 1 is no longer a suitable path of integration for ( p.2|) , because 
we have ReF(t) > for small t > 0, so that the endpoint t = is no longer dominant. 
What we need to do is to find an alternate contour of integration C running from t = to 
t = 1 such that ReF(t) < for all t G C except the endpoint t = 0. To see that we can 
indeed connect the endpoints t = and t = 1 by such a contour, let us look at the curve 
ReF(M,f) = in the complex t-plane. Its singular points are given by the simultaneous 
solutions of ReF(ti, f ) = and the set of equations 



dReF{u,v) 
du 

dReF{u,v) 
dv 



u-1 



= -y 



(u- l)2 + t;2 

V 







(B.15a) 
(B.15b) 



Now, the unique solution of (B.15) is given precisely by the saddle point (B.6); so singular 
points of the curve ReF(u,f) = can occur only when ReF(t^) = 0, i.e. when ^ E S (and 
in this case the unique singular point is the saddle point of F). Therefore, as we vary 
the curve Re F{u, f ) = can change its topological nature only when ^ crosses S. Now, no 
matter what the value of ^, the point t = 1 lies in the region Re F < 0, while the point t = 
lies on its boundary. And we have already shown that when ^ G ^1 U ^2, the points t = 
and t = 1 can be connected by a path (namely, the line segment < t < 1) lying entirely 
in the region ReF < except for the endpoint t = 0. And since we can get from A2 to 
anywhere in ^3 without crossing iS, it follows that for all ^ G ^3 there continues to exist a 
path C from t = to t = 1 lying entirely in the region ReF < except for the endpoint 
t = 0: see Figure P^a-c. Moreover, the same path C will work for some neighborhood of the 
initial ^. Therefore, the endpoint t = is again dominant, the standard theorems apply, and 
we obtain the same asymptotic expansion ([B.12| ); moreover, it is valid uniformly on compact 
subsets of A. 



We remark that in the smaUer domain |^| < 1, (B.12) can alternatively be proven by expanding the 
series representation 



(«)fc HY 



k=0 



ib)k k\ 



term-by-term in inverse powers of h. Slater |61, p. 66] mistakenly asserts that this holds for bounded ^, while 
Tricomi [|o[ pp. 111-112] correctly states that it works for j^j < 1. 
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Let us now look at the case when ^ E B. In this case, ReF(t^) > 0, so the saddle point 
dominates over both endpoints [namely, ReF(O) = and ReF(l) = —oo]. It suffices to find 
a path C running from t = through the saddle point and thence to t = 1, such that 
ReF(t) < ReF(ti,) for alH G C except t = t^. Making the change of variables t = ^{1 — t) 
[valid for ^ 7^ 0], we find that 

F{t) - F{Q = -{t- 1 -logt) , (B.16) 

so the curve ReF(t) = ReF{t^,) is precisely a dilated and rotated version of the Szego curve 
S. Moreover, the image of t = (namely, t = ^) lies in B, while the image of t = 1 (namely, 
t = 0) lies in Ai. So each of these two points can be connected to the saddle point t = t^. 



(whose image is t = 1) by a path whose image lies in B and Ai, respectively: see Figure p2 
This provides the needed path C. I 



B.2 Complex Zeros of iFi(a; b; b^) as 6 ^ 00 

We can now use some general results on the limit sets of complex zeros for certain 
sequences of analytic functions, proved recently by one of the authors [^], to deduce from 
the asymptotic expansions (p.l2| )/ ( |B.13|) the existence of zeros of iFi{a; b; bC,) near the Szego 
curve iS+ for large b. 

We begin with a "real" form of the theorem: 

Theorem B.2 Let D be a domain in C, and let G -D flR. Let (fm) be analytic functions 
on D, and let (am) be positive real constants such that d/ml""") uniformly bounded on 
compact subsets of D. Assume further that 

(a) For each x E D ClW, fm{x) is real and > 0. 

(b) For each x E D HM, lim am log fm{x) = g{x) exists and is finite. 

m— +00 

(c) g is not real-analytic at xq. 

Then, for all m sufficiently large, there exist zeros of fm such that lim z^ = Xq. 



Theorem B.2 is in fact well-known to workers in mathematical statistical mechanics: it is 



the contrapositive of an observation going back to Yang and Lee [^ concerning the genesis 



of phase transitions (see e.g. ||95|, Theorem 4.1, p. 51]). The proof of Theorem p.2| is, in fact, 
a direct transcription of the intuition sketched in Section ^. We first need a standard lemma 
from complex analysis: 

Lemma B.3 (log exp Vitali) Let D be a connected open subset of C", and let {gm}m=i 
be a family of analytic functions on D such that 

(a) For each compact C d D, sup sup Re (7^ (-2) < c)0 
and either 
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(bl) For some determining set S d D, lim gm{z) exists for all z & S 

m^oo 

or (h2) For some Zq E D and all multi-indices a, lim {D°'gm){zo) exists. 

Then gm converges, uniformly on compact subsets of D, to an analytic function g. 

We recall that S G D is called a determining set for D if g analytic on D and g \ S = 
implies g = 0. For example, any real environment (i.e. a set of the form zq + U, where U 
is a nonempty open subset of M") is a determining set. For n = 1, any set S having an 

P- 



accumulation point in Z) is a determining set. A proof of Lemma tB.3| can be found in 
343]; it is a special case of a vastly more powerful result, due to Montel (1912), on normal 
families of analytic functions |98| . 



Proof of Theorem |B . 2| . Suppose the contrary, i.e. suppose that there exists an e > 
and an infinite sequence mi < m2 < . . . such that none of the functions frm has a zero in 
the set = {z E C: \z — Xo\ < e} C D. Then, since is simply connected, log/^^ is 
analytic in (we take the branch that is real on fl M), and gi = (l/a^.) log/^^ satisfies 
limj^oo gi{x) = g{x) for x € fl M. Therefore, by Lemma [B.3| it follows that there exists an 
analytic function 'g on that coincides with g on D^d M. But this contradicts hypothesis 



The proof of Theorem |6.1| is now an immediate consequence of (6.13) [or ( [B.12| )/( [B.13| )] 
once we observe that |i-Fi(a; c; z)\ < e'^'^^l and that iFi(a; c; z) > for z real. 



Let us now state without proof a "complex" analogue of ( p.2| ) 



Theorem B.4 (Sokal [|93|| ) Let D be a domain in C, and let zq E D. Let (fm) be an- 
alytic functions on D, and let {am) be positive real constants such that (|/m|"™) o'^e uni- 
formly bounded on compact subsets of D. Suppose that there does not exist a neighborhood 
U 3 Zq and a function v on U that is either harmonic or else identically — oo such that 
liminf a^log |/m(z)| < v{z) < limsup log |/m(2;) | for all z G U. Then, for all m suffi- 

m ^oo m — >oo 

ciently large, there exist zeros of fm such that lim z^ = zq. 



The proof, which is not difficult, can be found in [93|; it uses elementary facts from the 



theory of normal families [^, ^ 



We can use Theorem [B.4| to deduce a vast extension of Theorem ^T] concerning the 
zeros of i-Fi(ci; 6; h^) for large h. Indeed, from the leading term of the asymptotic expansions 
(§;g)/(§^), we have 



lim -log|iFi(a;6;&OI 



for ^ G ^ 

Re(^-l-logO for^Gi3 



:B.17) 



We can then apply Theorem p.4| to any point ^ (which is the boundary between A 
and B) to deduce that: 

Theorem B.5 Fix a > and G 5+. Then there exist, for all sufficiently large b > 0, 
complex numbers zl satisfying i-Fi(a; h;zl) =0 and lim {zl/h) = C,o- 
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To the best of our knowledge, neither Theorem p3.5| nor its real specialization, Theo- 



rem 6J., can be found in the extensive literature on confluent hypergeometric functions 
(see H, [eg, HI H, [Tog, ig, H], IM], m, |10|] among many others). An old theorem of 
Tsvetkoff [^, Theorem 7] states that for b > 2a, every zero of i-Fi(a; 6; z) lies in the region 



Rez > 6 — 2a (see also Proposition 3.2] for a slight strengthening); Theorem |6]^ shows 
that this is asymptotically sharp as 6 ^ oo, to leading order. 



B.3 Asymptotic Expansion when ^ = 1 + 6 ^/^(^ 

In Section [B.l| we saw that the asymptotic expansion of i-Fi(a; b; b^) changes dramatically 
as ^ crosses the curve 5+, and in particular on the real axis as it crosses ^ = 1. We would 
now like to look more closely at the crossover behavior in a neighborhood of the point ^ = 1. 
It turns out that the right "magnification" is 6^/^: that is, we should look for an expansion 
of i-Fi(a; b; b^) when we take ^ = 1 + 6~^/^C with ( fixed and finite. 

Theorem B.6 Uniformly for ( in compact subsets of the complex plane, we have the asymp- 
totic expansion 



iFi (a; 6; 6(1 + 6-1/^0) 



m 



T{b-a) 

where Dy is the parabolic cylinder function 



Ai(a,C) A2(a,C) 
6V2 + b 



+ 0{b- 



-3/2N 



(B.18) 



:b.i9) 



and 



(B.20) 



a(a + l)Vt2) r(a + 4)/l a + 1 

I r(a + 6) 



/^-a-4(-C) 



(B.21) 



Corollary B.7 For large enough real b, iFi{a;b;b^) has complex zeros whose asymptotic 
expansion is given by 

^0 = 1 + ^ + ^ + ^ + 0(1/6^) (B.22) 
where Co is any zero of D^a{—0 , o-nd 

( , ^, /^-a-i(-Co) - |(a + 2)D_,_3(-Co) 
a = a{a + l) '\ (B.23) 
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+ 



a(a + l)2(a + 2) r(a + 4)/l a + 1 

7; ^-a-2{ — (,0> — 777— ^ (t + 



r(a) V4 



To prove Theorem [B.6| , we start with the integral representation 
^ = 1 + Expanding 



|ia2/)':,(-Co) + aaia + 1) [/^U_i(-Co) - |(a + 2)DV3(-Co)] 

^-a-4(-Co) 

(B.24) 
( |B.4| ) and set 

(B.25) 



I r(a + 6) 



= et + log(l-t) = (e-l)t - 5^ 



CO 



k=2 



k 



and making the change of variables s = b^^'^t, we can write ( p.2| ) in the form 
,F,ia;b;bil + b-'/\)) = 

{,1/2 



Cia,b)b'^/' J e-^/2+< (1-5^) 



-a-l 



^ exp 



CXD I. 

s 



E 

fc=3 



^5(fc-2)/2 



(B.26) 



The blow-up of the term (1 — s/6^/^)~"~^ as s t 6^/^ is harmless, as it is more than compen- 
sated (for large 6 > 0) by the suppression contained in the exp[— ■ ■ ■] factor [as can be seen 
by examining (p.l| ) as t t 1]. Therfeore, ( p.26|) becomes, for large b, 



iFi{a-b-b{l + b-^'\)) = C{a,b)b^'/^ J 6-''/^+''^ s^-Us [l + 0{b-^/^ 



C(a, b) 5-^/2 Y^a) e^'/^D_,(-C) [l + 0(6" 



1/2 



[B.27a) 



[B.27b) 



The corrections in powers of b can be easily extracted, and we arrive at expansion ( [B.18| ). 
Corollary p.7| then follows using the implicit function theorem. | 



C Correlation Inequalities 

In this appendix we prove some correlation inequalities for A^-component a-models with 
Hamiltonian of the form 

^(M) = -$^v.,K.-a,), (c.i) 

where each spin cr^. is a unit vector in and each Vxy is a real- valued function on [—1,1]. 
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C.l Tools #1: Griffiths Inequalities for the Ising Model 

Consider an Ising model consisting of finitely many spins Sx = ±1, with Hamiltonian 



H{{e}) = - E ^^^^ • 

A 

Here we have used the notation = H^x""' where A = {A^} is an arbitrary multi-index 

X 

(i.e. a vector of non- negative integers). Of course, for the Ising model it suffices to use 
Ax = 0, 1; the sum in (|C.2| ) runs over all such multi-indices. Thus, we have allowed quite 
general multi-spin interactions. Of course, magnetic fields and pair interactions are included 
as special cases. 

We then have the following well-known inequalities: 

Theorem C.l (Griffiths' first and second inequalities) Assume that Ja > for all 
A. Then: 

[G-I] (e^) > for all A. 

[G-II] {e^;s^) = {e^e^) - {e^){e^) > for all A, B. 

Theorem C.2 (Griffiths' comparison inequality) Assume that J a > \J'a\ for all A. 
Then 

[G-Comp] > for all A. 

Elementary proofs of (vastly generalized versions of) Theorems |C.1| and |C.2| can be found in 



|103|| and [|104|| , respectively. 



C.2 Tools #2: Ginibre Inequalities for the XY Model 

Next consider a classical XY (plane-rotator) model consisting of finitely many spins 
9x € [0,27r), with Hamiltonian 

H{{e]) = -^Re{JAe'^') (C.3a) 

A 

= - ^ |Ja| cos(v4 ■ 6* + arg J^) , (C.3b) 

A 

where the sum runs over all vectors A = {A^} of (signed) integers, the interactions J a are 
arbitrary complex numbers, and A ■ 6 = A^Ox- We then have the following inequalities: 

Theorem C.3 (Ginibre's first and second inequalities) Assume that Ja > for all 
A. Then: 

[Ginibre-I] {cos A -9) > for all A. 
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[Ginibre— II] (cos A-9; cos B- 
for all A,B. 



{cos A-9 cos B -9) - {cos A- 9) {cos B ■ 9) > 



Theorem C.4 (Ginibre comparison inequality) Assume that > {J'y^l for all A. Then 

[Ginibre— Comp] {cosA-9)j > \{cos{A ■ 9 + (fA)) j'\ for all A and all (real) 
angles (f^. 



Proofs of (generalized versions of) Theorems p.3| and |C.4| can be found ||105|| and ||106 
respectively. 



C.3 Correlation Inequalities for ^ cr-Models 

We now prove some correlation inequalities for A^-component cr-models. Our method — 
which is in essence due to Frohlich ||107|| — is to embed an Ising or XY model into the 
A^-component model, and then apply the Griffiths or Ginibre comparison inequality to this 
embedded model (which has annealed random couplings) 

We consider A^-component cx-models of the type 

H{{ct}) = - ^V.,(a.-a,) (C.4) 

(and even more general models). It is convenient to split Vxy into its even and odd parts: 

-s) 



xy\ 



(C.5) 



We begin by proving that correlations in the cr-model are bounded above by the correspond- 
ing correlations in the Ising model: 

Theorem C.5 Consider a a-model defined by the Boltzmann-Gibbs measure 
Z~^e^^^^^^^dfj,{{a}), where H is defined by ( \C.4i ), and c//i({cr}) is an arbitrary nonneg- 
ative measure that is invariant under the gauge transformations — > Tj^cTx (Vx = ^^)- 
Then, for any set of sites xi, yi, . . . , Xn, y-n, we have 



""j/n ) ) -H",/^ I — {^x\^y\ ' ' ' ^Xn^yn) Ising,J {\^ x\' 



yi I 



— {^xi^yi ' ' ' ^Xn^yn) Ising,J ((""i-'i ' 



yiy 



^ {^xi^yi ■ ■ ' ^Xn^yn) Ising,J j 



where 



Jxy 



IV, 



xy\ 



sup |V_(s)| . 



se[-i,i] 



'x„ - (^y„\)H,ti 

(C.6a) 

{(^x„- (TyJ )H,t, 

(C.6b) 
(C.6c) 

(C.7) 



•^^ A slightly weaker version of this method was applied earlier by Mack and Petkova |10S|. They were 
apparently unaware of the Griffiths and Ginibre comparison inequalities, so they relied instead on the 
(somewhat weaker) Griffiths II and Ginibre II inequalities. To satisfy the ferromagnetism hypotheses of 
these latter inequalities, they had to impose a "no- vortex" (or "no-monopole" ) constraint. This constraint 
is now seen to be optional but not necessary. 
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Proof. The proof is based on the trivial identity 

I F{{cT})df^i{cT}) = F{{ecT})df^{{cT})dui{e}), (C.8) 

where di'{{e}) is normahzed counting measure on the configurations of Ising spins = ±1, 
and F is an arbitrary function. This identity expresses the Z2 gauge invariance of dfi{{(T}). 

Without loss of generality we can assume that each V^y is an odd function, since the even 
parts of Vxy can be absorbed into /z. 



Let us apply ( |C.8| ) to F{{(t}) = (cr^.^- a-yj ■ ■ 



-H{{(T}). 



^xi^yi ' ' ' ^Xn^yn \^xi ' ^yij 



rf/i({cr}) (cr^.,- cr 



[(^xr.-(TyJ I ex,ey,---e,„ey^ Y[e'^'y^^y^'^^-'^yUu{{e}) 



(xy) 



(C.9) 

We now perform the integration over {e}; by the Griffiths comparison inequality (Theorem 
U.2[) we have 



^xi^yi ■ ■ ■ ^Xn^ 



n 

(xy) 



du{{e}) 



< {e.,ey,---e,,^ey„)j..ny,j j lle^-^y''^y^'^-'^yUu{{e}) 

{xy) 



(C.IO) 



It follows that 



< {e^.ey, ■ ■■e^^EyJlsingJ / dfi{{(7}) \ {(Tx,-(7 



yiy 



[o- xn ■ cr y„ 



(xy) 



\^xi^y\ 



' ^x„£y„) 



Ising, J 



['^xi ' '^yij 



(C.ll) 



The rest follows from the Schwarz inequality together with the trivial bound |cr| < 1. 



Corollary C.6 For the mixed isovector /isotensor model 

K 



H 



E 

(xy) 



Jxy^x ' 



<^x - <^y) 



we have 



\{{(Tx^- (Ty^) ■ ■ ■ {(Tx„- CryJ)j^K\ < {^xi^yi- ■ ■£x„£yjlsing,\j\ ■ 



(C.12) 



(C.13) 
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Remarks. 1. The bound of Theorem |C.5| and Corollary |C.6| holds for arbitrary products 



but it is only of interest for isovector correlation functions. For isotensor 



correlation functions (i.e. those in which each spin cr^ appears an even number of times), 
the corresponding Ising- model correlation is 1, and the bound is trivial. 

2. The only restriction on the measure n is that it be Z2-gauge- invariant. An interesting 
family of examples is given by 



dfi{{(T}) = exp 



/^pj]^ n sgn(cr^-o-y) 

5 (xy)GS 



(C.14) 



where the sum runs over plaquettes S. For l3p = 0, this is the ordinary cr-model. For /3p — > 
+00, this is the no-vortex cx-model ||108| , |109|| , in which the constraint n(a;y)G5 ^x- (^y ^ ^ 
imposed for every plaquette S. 

3. This proof is based on the embedding {e,cr) i— ecr. More general embeddings could 
be considered, i.e. {e, cr) {I — P)(t + sPcr where P is any orthogonal projection in M^. 
(If rankP = 1, this is the Wolff embedding if rankP = k, this is the co dimension- /c 

embedding discussed in |^T|]; if rankP = A^, it is the embedding used above.) We do not 
know whether embeddings with rank P < N might be useful for proving correlation inequal- 
ities. We note, however, that the Wolff embedding has been employed for analytic purposes 
(rigorous lower bounds on the correlation function via a Fortuin-Kasteleyn representation) 



by Patrascioiu and Seller |110| and Aizenman [111]. 



Next we consider bounds relating the A^-component a-model to an XY model. For 
simplicity we consider first an ordinary A^-vector model 



and compare it to an XY model 



i^y) 



xy 



COS 



y ' 



(C.15) 



(C.16) 



Theorem C.7 Consider an N-vector model defined by the Boltzmann-Gihhs measure 
Z^^e^^^^^^^diJL{{(T}), where H is defined by ( \C.13i ) and dn{{(T}) = Yl di^{cr r,.) . Consider 

X 

for comparison an XY model defined by the Boltzmann-Cibbs measure Z^^ e^^'^^^^'^ dv{{6}) , 
where H is defined by ( \C.1(^ and c/z/({cr}) = J^fi^a;- Assume that \ J'^y\ < Jxy for all x,y. 

X 

Then 

\{crx- (Ty)N,j'\ < {cos{9x - 9y)) xY,j (C.17) 



and 



1 



1 



{N-l){N + 2) 



N [2N -{N- 2)(cos2(^, - ^,))xy,j] 
(cos2(6'a. - 9y))xY,j ■ 



1 

N 



(cos 2(6*^ - Oy))xY,j 
(C.18a) 
(C.lSb) 
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Proof. Let R{6) e SO{N) be the rotation through an angle 9 in the 1-2 plane. We will 
employ an XY embedding {6^, cr^) i-^ R{9x)(Tx. The proof is based on the trivial identity 

j F{{a})di^{{<T}) = F{{R{e)cT})dfi{{a})dum) , (C.19) 

where dfi{{(T}) is uniform measure on the product of spheres, di'{{6}) is uniform measure 
on the product of circles^ and F is an arbitrary function. 

Define al = and = . . . , ^). Let us apply ([UJQD to F({cr}) = 

(o-|l-a;)e~^«'^»: 



cr" ■ cr 



J)e-^«'^»dM{^}) 



J J cos% - 9y + Z(^i„ ^5)] (1 - - aiY' e-^«^(^)'^» df^{{a}) dum) 



(C.20) 



where Z{(Tx, cr|) denotes the angle between cr| and CTy. Note that if ({-R(6')cr}), considered as 

a function of {^} for fixed {cr}, is an XY model with couplings J^y = Jxy\f^x\ \cr\\e^^''^'^'^y^ . 

Since \Jxy\ < \J'xy\ — Jxyi "wc cau perform the integration over {^}, using the Ginibre 
comparison inequality (Theorem |C.4|) to obtain an upper bound: 



|cos[0.,-^^, + Z(a|.,a;)]e-^(™'^» 



dum) 



< (cosf 



(C.2i: 



It follows that 



(cos(^. - ^,))xy,. J )^/^(l-<Y/^e-^«'^»d/i({'^}) • 



(C.22) 



In terms of normalized quantities, this says that 

\{crl^l)N,j'\ < {cos{e.^ -ey))xY, J {{I -^iY'il-^fy^')^^, ■ (C.23) 

Now we use the inequality ab < |(a^ + b^) together with global 0{N) invariance to conclude 
that 



2 

iV 



(C.24) 
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Hence 



{(tI- (tI)n,j'\ < j;^ {cos{e^-ey))xY,j ■ 



(C.25) 



Using 0{N) invariance again in the form (cTx- (Ty) = {N/2){(t\ - cr|), we complete the proof 

of (013) • 

Next let us apply (|C.19|) to F({cr}) = [cri - (T\Ye^^''^'-''^\ By the same steps as before, 
we show that 



Now we exploit 0{N) invariance in the form of identities 

(Ti-^)) = 



N,J' 



(C.26) 

(C.27a) 
(C.27b) 



where the isotensor field tJ"^^ is defined in ( p.7| ). After some algebra we obtain 



N,J' 



so that (|C26D can be written as 



+ 



QN -A 



fq(^fq2 J^N -2) 
N-2 



GT{x,y) 



N 



N-2 



N 



+ 



AN -8 



N{m + N -2) 



GT{x,y) 



(C.28a) 
(C.28b) 
(C.28c) 



+ 



QN -A 



m Nim + N -2) 



+ 



AN -8 



m N{m + N -2) 



GT{x,y) 



(C.29) 



where ■jxy = (cos2(^^a; — 6y))xY,.j- The bounds (CIS) follow by simple algebra. 



Remarks. 1. Note that the additive term in (C.18) is just what is needed to give 
the (properly truncated) isotensor correlation function, i.e. 



Gt{x, y)N,j' < 1 



N 



(cos 2(6*^ - Oy))xY, 



J ■ 



(C.30) 



2. A partially alternate proof of Theorem |C.7| can be based on the parametrization 
cTj. = ^(1 — cr^^y^'^ cos9x, (1 — (T^^y^^ sin 6*2,, cr^^. Thus, 9x is here the polar coordinate of 

crl (whereas it was previously an additional rotation applied to cr|.). Under the identification 



the uniform measure d^{{a}) becomes di'{{6}) dfi{{a-^}) for a suitable 
measure Jl (whose form is of no interest to us). Now the Hamiltonian H{{cr}), considered 
as a function of {9} with {cr-*-} fixed, has couplings Jxy = ^L(l — cri"^)^^^(l — (^x'^Y^'^i ^^'^ 



56 



so is ferromagnetic if the original model is (i.e. if J'^y > 0). The remainder of the proof is as 
before, but in the ferromagnetic case we can now use the Ginibre II inequality rather than 
the stronger Ginibre comparison inequality. However, the first proof has the advantage that 
it generalizes readily to any cr-model whose invariance group contains a U{1) subgroup, e.g. 
to principal chiral models. 

2. One could consider different embeddings of XY spins into A^-vector spins, e.g. different 
one-parameter subgroups R{0) G SO{N). We do not know whether such embeddings might 
be useful for proving correlation inequalities. 

Finally, let us consider the iV-component mixed isovector/isotensor model 

H{{^}) = - E [J'^y^^-^y + iKyi^^-^y?] , (C.31) 

and compare it to a generalized XY model 

Hm) = -Yl [^-j/ - + ^-y - ey)] . (c.32) 



Theorem C.8 Consider an N -component mixed isovector/isotensor model defined by the 
Boltzmann-Gibbs measure Z^^e^^^^^^^djj.{{(T}) , where H is defined by ( \C. 31 ), and dn({(T}) = 
Yl dQ{(Tx)- Consider for comparison a generalized XY model defined by the Boltzmann-Gibbs 

X 

measure Z^^e^^'^^^^'>du{{6}), where H is defined by ( \C. 32J , and du{{(T}) = Y[dOx- Assume 

X 

that 

I /' I if \K' \ < \ I' 

I XV I V \-^^ XV I — \ XI 



xy\ J \ xy\ — \" xy\ 

\(\ + \Iky[S\' jf\K' I > I 7' I 
4K^xyl ^ \K'^y\) y \^^xy\ ^ \^xy\ 



< Jxy (C.33) 



and 



\\Ky\ < K.y (C.34) 
for allx,y. Then inequalities ( \C. 1 ?| j and (C.18) hold for {■)n,j',k' o^nd {■)xy,j,k- 

Remark. The left-hand side of ( |G.33| ) is bounded above by | J^^^l + So one can 

get a simpler (but weaker) version of the theorem by using this formula in the hypothesis. 

Proof. The proof is identical to that of Theorem |C.7| . We find that H{{R{6)cr}), considered 
as a function of {6} for fixed {cr}, is a generalized XY model with couplings 

ly = (J:, + ir>^0|all||a;|e^^('^-'^il) (C.35a) 

Kxy = ^i^;,kErk!HV^^('^^''^^) (C.35b) 

By straightforward calculus we see that \Jxy\ < the left-hand side of (|C.33|) . So we can use 
Ginibre comparison inequality (Theorem |U.4| ) as before. The rest of the proof concerns the 
treatment of the various observables (as opposed to Hamiltonians) , and so is unchanged. ■ 
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Figure 1: Phase diagram of the RP^~^ model at = oo in the (/3t, (i)-plane. Transition 
from white-noise phase to isotropic phase is second-order for < d < 3/2 and first-order 
for 3/2 < d < dt, ^ 2.38403. Transition from isotropic phase to nematic ordered phase is 
second-order for 2 < c? < d*. Transition from white- noise phase to nematic ordered phase is 
first-order for d> d*. 
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Figure 2: .^0(7) versus 7 for various values of Pt, for (a) d — 1, (h) d — 3/2, (c) d — 2. 
Since .^^0(7) = we show only 7 > 0. 
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Figure 3: .^0(7) (solid curve) and To*{f) (dashed curve) plotted versus / for various values 
of Pt: for (a) d — 9/4, (b) d = d^, ^ 2.38403, (c) d — 3. The junction between the solid and 
dashed curves corresponds to 7 = l/d, / = 
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Figure 4: d'^c plotted versus d for 3/2 < d < d^, 2.38403. Here 7c is the 7 value at the 
first-order phase transition of the RP°° model [see (|4.6| )]; it tends to zero as (i | 3/2 and to 
1/ d a.s d } dt,. 




Figure 5: JF^ = — (/* — l)/4 + (log/* + fi'*)/2 plotted versus d. We have JF,, = at the 
dimension ^ 2.38403. 
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Figure 6: Pt versus 7 at /3y = 0, for (a) d = 1, (b) d = 3/2,Jc) d = 2, (d) d = 9/4, 
(e) d = 2.38403, (f) d = 3. For 3/2 < d < J fa 2.55031, Pt reaches a minimum at 
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Figure 7: Free energy T versus fix at fiy = 0, for (a) d = 1, (h) d = 3/2, (c) d = 2, (d) 
d = 9/4, (e) (i = ~ 2.38403, (f) d = 3. In (d)-(f), the dashed curve corresponds to the 
condensate phase. 
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Figure 8: Phase diagram in the /3t) -plane for dimensions < d < 3/2. There is a 
second-order transition at {(3v,Pt) — (0, 1) and a hne of first-order transitions at (3v — 0, 
Pt>1. 




Figure 9: Phase diagram in (/^y, /9T)-plane for dimensions 3/2 < d <2, shown here for the 
case d = 2. The first-order transition curves (3t,c{(3v) are almost straight fines but are in fact 
sfigfitly convex. In d = 2, the triple point lies at {Pv, Pt) ~ (0, 0.956) and 7c ~ 0.487, while 
the endpoints of the first-order transition curves lie at {±Pv,e: PT,e) ~ (±0.0134, 0.905) and 
7e 0.411. 
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Figure 10: Phase diagram in the (/Sy, /9t) -plane for dimensions 2 < < rf* ^ 2.38403. 
In addition to the first-order transitions present for 3/2 < < 2 at (3t,c{13v)i there is an 
additional pair of second-order transition hnes (dashed) at Pt,c'{Pv) given by (|5.1|). The 
phase at the upper right (resp. upper left) of the diagram exhibits ferromagnetic (resp. 
antiferromagnetic) long-range order. 
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Figure 11: Phase diagram in the (/^y, /3t) -plane for dimension d = ^ 2.38403. The 
two transitions coincide at Pv = 0. The phase at the upper right (resp. upper left) of the 
diagram exhibits ferromagnetic (resp. antiferromagnetic) long-range order. 
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Figure 12: Phase diagram in the /Sy) -plane for dimensions ^ 2.38403 < d < d„ = 
3. Second-order transitions (dashed hnes) occur only where they lie above the first-order 
transition curves; this occurs for |/5y| > Pv,cc'- The phase at the upper right (resp. upper 
left) of the diagram exhibits ferromagnetic (resp. antiferromagnetic) long-range order. 
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Figure 13: Phase diagram in (/3y, /9T)-plane for d > d** = 3, shown here for the case 
d — 3. The first-order transition curves Pt,c{Pv) are almost straight lines but are in fact 
slightly convex. In = 3, the triple point lies at (/Sy, /^t) ~ (0, 0.73487), while the endpoints 
of the first-order transition curves lie at {f3v,e, Pr.e) = (='=/*(3)[2 — /^,(3)]/6, /*(3)^/6) ~ 
(±0.122224, 0.383238). The phase at the upper right (resp. upper left) of the diagram exhibits 
ferromagnetic (resp. antiferromagnetic) long-range order. 
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Figure 14: Coordinate /3T,e of the endpoint of the first-order transition curve, as function 
of dimension d. For d > c?** = 3 it is given by (3T,e = fl/'^d (dashed fine) and tends to zero 
as d — > oo. 




Figure 15: Coordinate I3v,e of the endpoint of the first-order transition curve, as function 
of dimension d. For d > d^^ — ii is given by (3v,e — Pv,2 = /*(2 — f*)/'^d (dashed line) and 
tends to zero as d — > oo. 
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Figure 16: "Notations guide" for the parametric plots of JF versus (3t- Points 1 and 2 
are, respectively, the local maximum and minimum of 0t versus /. Point c is the first- 
order transition at Pt,c{Pv)- Point c' is the second-order normal-to-condensate transition at 
f — f*. The condensate phase / > /* is shown as a dashed line. 
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Figure 17: (3t versus /, and JF versus (3t-, for d = 2 at (a) (3y = 0, (b) jSy = 0.001, 
(c) = 0.005,^(d) 3y = ^v,e - 0.0134 [the dot indicates the endpoint (py^e^T^e)], (e) 
= 0.017, (f) = 0.02. 
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Figure 18: (5t versus /, and T versus Pt, for d = 2.25 at (a) j3y = 0, (b) j3y = 0.01, 
(c) = 0.02, (d) py = 0.032, (e) py = py^e - 0.03464 (note that pT,e - 0.795978 and 
7e ^ 0.416596), (f) I3y = 0.035. 
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Figure 19: Pt versus /, and versus Pt, for d = 2.45 at (a) Py = 0, (b) Py = 0.01, (c) Py = 
py^^^, ^ 0.0237503, (d) py = 0.04, (e) py = py^^ ~ 0.057583 (note that pT,e ~ 0.69010393 
and 7e ^ 0.399979), (f) py = 0.07. 
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Figure 20: Pt versus /, and T versus /3t, for d — 2.75 at (a) Pv — 0, (b) Pv — 0.05, (c) 

py = py^^ ^ 0.0892914 (d) Pv = Pv^cd ^J).0946534, (e) Pv = Pv,e - 0.0965845 (note that 
pT^e - 0.5158930 and 7e ^ 0.363383), (f) pv = 0.1. 
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Figure 21: Blow-ups of the parametric plots for d — 2.75. (a) At /3y = /3v^,2 ~ 0.08929, 
points c' and 2 coincide, (b) At — ^v,cc' ~ 0.09465, points c' and c coincide, (c) At 
/3y = /3v^,e ~ 0.09658, the "triangle" disappears; this occurs within the normal phase. 
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Figure 22: (3t versus /, and versus Pt, for = 3 at (a) Pv = 0, (b) /3y = 0.025, 
Pv = 0.05, (d) pv = Pv,2 = Pv,cc' = Pv,e ~ 0.122224, (e) pv = 0.15. 
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Figure 23: Pt versus /, and JF versus for d = 5 at (a) Pv = 0, (b) Pv = 0.01, 
py = py^^ K 0.0571384, (d) pv = 0.075, (e) pv = Pv,2 = Pv,cc' = Pv,e ~ 0.0975568, 
Pv = 0.1. 
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Figure 25: /3v^,e) /3y,cc'; /3y,2 and /3v^,i as functions of dimension d. The curve /3y,e smoothly 
reaches at d = 3/2; the curve /3y,cc' hits at d = 2.38403; the curve /3y_2 = /*(2 — 
/*)/(2d) hits at (i = J ^ 2.55391; the curves /9y,e, /9v",cc' and /9y,2 all merge at d = d** = 3; 
the curve fiv,\ branches off the /3y,e — l^v,cc' — Pv,2 curve at d — 4, and then approaches zero 
as c? — >• oo. 
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Figure 26: Approach to the N = oo hmit in the one-dimensional RP^ ^ model: curves 
shown are = 10, 25, 100, 1000, oo. (a) Free energy F{(3) = -(1/A^) logZjv(iV/3). (b) 
Correlation decay rate vn^2 = e"™''^ where rriT is the isotensor mass gap. 
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Figure 28: Magnitude of the first imaginary zero of Zj^{N jSy^N Pt) plotted versus A^~^/^, 
for (3t = 0, 1/4, 1/2, 1. Solid box indicates the exact limiting value ( [7. 221 ) for = oo. Slope 
of straight hne is given by a best fit to the data for > 50, except for /?t = 0, when it is 
given by the exact asymptotic formula (see footnote p^). 
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Figure 29: Conjectured phase diagram in the /3T)-plane for the mixed isovec- 
tor/isotensor model with 2 < A?" < oo in dimension d > 2. First-order transitions are 
shown as sohd hnes, and second-order transitions as dashed hnes. 
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Figure 30: Szego curve x"^ + y'^ = e^^^~^^ in the complex ^ — x + iy plane: we denote by (S_ 
(resp. 5+) the part of this curve at x < 1 (resp. x > 1). The regions Ai, A2, A3 and B are 
indicated. 
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Figure 31: Curve ReF(i) = in the complex t — u^iv plane, for (a) ^ G Ai-, (b) ^ on the 
boundary between A2 and ^3, (c) ^ e ^3, (d) ^ e «S+, (e) ^ G S. The points t — ^ and i = 1 
can be connected by a path lying in the region Re F < whenever t ^ but not otherwise. 
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Figure 32: Curve ReF(t) — ReF{t^) in the complex t — u + iv plane, for a point ^ G B. 
One possible path C is shown in dashes: it runs from t = to the saddle point t^, and thence 
to t = 1, staying entirely in the region ReF{t) < ReF{t^) except at t = t^. 
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Figure 33: Szcgo curve x'^ + y"^ — e^^^ together with some zeros of iFi(a; b; b{x + iy)) for 
a — 1/2, h — 1000. Note the accumulation of zeros densely on 5+. 
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